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SUMMARY 


This report documents the derivation and definition of a linear aircraft model for a rigid aircraft of constant 
mass flying over a flat, nonrotating earth. The derivation makes no assumptions of reference trajectory or 
vehicle symmetry. The linear system equations are derived and evaluated along a general trajectory and 
include both aircraft dynamics and observation variables. 


INTRODUCTION 


The need for linear models of aircraft for the analysis of vehicle dynamics and control law design is well 
known. These models are widely used, not only for computer applications but also for quick approximations 
and desk calculations. Whereas the use of these models is well understood and well documented, their 
derivation is not. The lack of documentation and, occasionally, understanding of the derivation of linear 
models is a hindrance to communication, training, and application. 


This report details the development of the linear model of a rigid aircraft of constant mass, flying over a 
flat, nonrotating earth. This model consists of a state equation and an observation (or measurement) equa- 
tion. The system equations have been broadly formulated to accommodate a wide variety of applications. 
The linear state equation is derived from the nonlinear six-degree-of-freedom equations of motion. The 
linear observation equation is derived from a collection of nonlinear equations representing state variables, 
time derivatives of state variables, control inputs, and flightpath, air data, and other parameters. The linear 
model is developed about a nominal trajectory that is general. 


Whereas it is common to assume symmetric aerodynamics and mass distribution, or a straight and level 
trajectory, or both (Clancy, 1975; Dommasch and others, 1967; Etkin, 1972; McRuer and others, 1973; 
Northrop Aircraft, 1952; Thelander, 1965), these assumptions limit the generality of the linear model. The 
principal contribution of this report is a solution of the general problem of deriving a linear model of a rigid 
aircraft without making these simplifying assumptions. By defining the initial conditions (of the nominal 
trajectory) for straight and level flight and setting the asymmetric aerodynamic and inertia terms to zero, 
one can easily obtain the more traditional linear models from the linear model derived in this report. 


Another significant contribution of this report is the derivation and definition of a linear observation 
(measurement) model. The observation model is often entirely neglected in standard texts. A thorough 
treatment of common aircraft measurements is presented by Gainer and Hoffman (1972), and Gracey (1980) 
provides a detailed discussion of speed and altitude measurements. However, neither of these references 
present linear models of these measurements. This report relies heavily on these two references and uses their 
results as one of the bases for the nonlinear measurement equations from which the linear measurement 
model is derived. Also included in this report is a large number of other measurements or variables for 
observation that have been found to be useful in vehicle analysis and control law design. 


Duke and others (1987) describe a FORTRAN program called LINEAR that derives a linear aircraft 
model by numerical differencing (Dieudonne, 1978). The program LINEAR produces a linear aircraft model 
(both state and observation matrices) that is equivalent to the linear models defined in this report, 


This report is divided into two main sections that define the reference systems and nonlinear state and 
observation equations (section 1) and derive a linear model presented in the appendixes (section 2). The 
appendixes contain a definition of the linear aerodynamic model used in this report (app. A), a derivation 
of the wind axis translational acceleration parameters (app. B), generalized linear derivatives of the non- 
linear state and observation equations (app. C), and the individual derivatives of the state and observation 
equations (app. D). The details of the principal results of this report are presented in appendix D. 
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total aerodynamic axial force, lb 

speed of sound, ft/sec 

normal accelerometer output, g¢ 

output of normal accelerometer not at vehicle center of gravity, g 

output of accelerometer aligned with vehicle body z axis, g 

output of accelerometer aligned with body z axis, not at vehicle center of gravity, g 
kinematic acceleration in vehicle body z axis, g 


output of accelerometer aligned with vehicle body y axis, g¢ 

output of accelerometer aligned with body y axis, not at vehicle center of gravity, g 
kinematic acceleration in the vehicle body y axis, g 

output of accelerometer aligned with vehicle body z axis, g 

output of accelerometer aligned with body z axis, not at vehicle center of gravity, g 
kinematic acceleration in vehicle body z axis, g 

reference span, ft 

generalized force or moment coefficient 

derivative of generalized force or moment coefficient with respect to arbitrary variable z 
reference aerodynamic chord, ft 

total aerodynamic drag, Ib 
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specific energy, ft 

arbitrary force or moment 

flightpath acceleration, g 

acceleration due to gravity, ft/sec? 

acceleration due to gravity at sea level, ft/sec? 
altitude, ft 

altitude measurement not at vehicle center of gravity, ft 
inertia tensor 

moment of inertia about z body axis, slug-ft? 
product of inertia in x-y body axis plane, slug-ft? 
product of inertia in z—z body axis plane, slug-ft? 
moment of inertia about y body axis, slug-ft? 
product of inertia in y-z body axis plane, slug-ft? 
moment of inertia about z body axis, slug-ft? 
II, — I2, 

Tayls + Lyzlez 

Ieylyz + Lylrz 


Inlz — I2, 
Telyz + Leylez 
Inly — H3, 


total moment about x body axis, ft-lb; or, total aerodynamic lift, Ib 
unit length, ft 
total moment about y body axis, ft-lb; or, Mach number 


DWP AA 


‘S. 
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vehicle mass, slugs 


total moment about z body axis, ft-lb; or, total aerodynamic normal force, |b 


load factor 

specific power, ft/sec 

roll rate (about x body axis), rad/sec 
static or free-stream pressure, lb/ft? 
stability axis roll rate, rad/sec 

total pressure, lb/ft? 

pitch rate (about y body axis), rad/sec 
dynamic pressure, lb/ft? 

impact pressure, lb/ft? 

Mach meter calibration ratio 

stability axis pitch rate, rad/sec 
Reynolds number 

Reynolds number per unit length, ft-? 
yaw rate (about z body axis), rad/sec 
stability axis yaw rate, rad/sec 
surface area of wing, ft? 


total angular momentum; or, ambient or free-stream temperature, oR 


total temperature, °R 

time 

velocity along x body axis, ft/sec 

vehicle velocity, ft/sec 

velocity along y body axis, ft/sec 

velocity along z body axis, ft/sec 

total aerodynamic force along z body axis, lb 
total gravitational force along x body axis, lb 
total thrust force along x body axis, Ib 
vehicle position along x earth axis, ft 

total aerodynamic sideforce, lb 

total aerodynamic force along y body axis, lb 
total gravitational force along y body axis, |b 
total thrust force along y body axis, lb 
vehicle position along y earth axis, ft 

total aerodynamic force along z body axis, lb 
total gravitational force along z body axis, lb 
total thrust force along z body axis, lb 
vehicle position along z earth axis, ft 

angle of attack, rad 

angle-of-attack measurement not at vehicle center of gravity, rad 
angle of sideslip, rad 

angle-of-sideslip measurement not at vehicle center of gravity, rad 
flightpath angle, rad 

ith control surface deflection 

pitch angle, rad 

coefficient of viscosity, lb/ft-sec 

density of air, lb/ft? 

arbitrary function 


Vectors 
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Matrices 


bank angle, rad 
heading angle, rad 


body axis acceleration vector 

attitude vector of Euler angles 

total force vector 

state vector function 

observation vector function 

total angular momentum vector 

sum of higher order terms in Taylor series 
total moment vector 

position vector in earth axis system 

input or control vector 

vehicle velocity vector 

state vector 

observation vector 

perturbation of control vector 
perturbation of state vector 

perturbation of time derivative of state vector 
rotational velocity vector 


state matrix of the generalized state equation, Cx = Ax+ Bu 

state matrix of the state equation, x = A’x + B’u 

control matrix of the generalized state equation, Ck = Ax + Bu 

control matrix of the state equation, x = A’x + B’u 

system matrix of the generalized state equation, Ck = Ax+ Bu 

feedforward matrix of the generalized observation equation, y= Hx+Gx+ Fu 
feedforward matrix of the observation equation, y = A’x+ Fu 

derivative observation matrix of the generalized observation equation, y= Hx+Gx+ Fu 
observation matrix of the generalized observation equation, y = Hx+Gx+ Fu 
observation matrix of the observation equation, y = H’x+ F’u 

intertia tensor 

scaling matrix for inertia tensor 

transformation matrix from earth to body axes 

transformation matrix from earth to body axes 

angular velocity matrix in the generalized state equation, Tx = f[x(t),x(t), u(t)] 
n Xm matrix of 0 values 

an n X m matrix with values of 1 on the diagonal 


Subscripts 


a aerodynamic; or static or, free stream 
b body axis system 

D drag 

g gravitational 

h displacement of altitude instrument 

h displacement of altitude rate instrument 
yt not at vehicle center of gravity 
kinematic 

lift 

rolling moment 

pitching moment 

yawing moment 

orthogonal 

power plant induced 

stability axis; or, specific 

thrust 

total 

vehicle-carried vertical axis system 
wind reference axis system 
displacement in x body axis 

z~-y body axis plane 

x—z body axis plane 

sideforce 

displacement in y body axis 

yz y-z body axis plane 

z displacement in the z body axis 

0 at sea level, standard day conditions; or, nominal conditions 


= 
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eme ge ect yes gee 


Superscript 


T transpose 


1 NONLINEAR SYSTEM EQUATIONS 


The motion of an aircraft as a rigid body can be described by a set of six nonlinear simultaneous second- 
order differential equations. These equations, representing the translational and rotational motion of the 
vehicle, can be formulated in the notation of Kwakernaak and Sivan (1972) and Dieudonne (1978) as a 
time-invariant system expressed as 

x(t) = f[x(4), u(¢)] (1-1) 
where x(t) is the 12-dimensional time-varying state vector (¢ being time), x(t) is the derivative of x(t) with 
respect to time, u(¢) is the k-dimensional time-varying input or control vector, and f is a 12-dimensional 
nonlinear function expressing the six-degree-of-freedom rigid body equations. 


Measurements of the vehicle state can be represented by the observation equation 


y(t) = g[x(z), u(z)] (1-2) 


where y(t) is an ¢-dimensional time-varying observation vector and g is an ¢-dimensional nonlinear func- 
tion expressing the relationship of the true vehicle state and control vectors to the observed parameters. 
Typically, the function g characterizes the dynamics and location of the sensors. 


For the aircraft analysis and design problem, both the nonlinear and linear system equations are formu- 
lated more broadly than just described (Edwards, 1976; Maine and Iliff, 1980, 1986). The nonlinear system 
equations include x(t) terms in both the state and observation functions. In fact, in the most extended 
form the state equation is expressed in terms of transformed variables (discussed in section 1.2.1). These 
generalized equations form the basis of the analysis in this report. The generalized system equations are 


T(t) = f[x(t), X(t), u(2)] (1-3) 


y(t) = glx(t), x(£), u(t) (1-4) 


where T is a constant 12 x 12 angular velocity matrix. 


1.1 Definition of Reference Systems 


While numerous reference systems are used in aerospace applications, this report is limited to four reference 
systems: the body, the wind, the vehicle-carried vertical, and the topodetic reference systems. The stability 
axes are also defined even though this reference system is used only to define the stability axis rotational 
rates (section 1.3.8). 


Within this report the translational equations are referenced to the wind axes, and the rotational 
equations are referenced to the body axes. Measurement equations are primarily referenced to the body 
axes when the use of a reference system is needed. The use of this mixed axis system definition in both 
the nonlinear and linear models is related to the measurability and meaningfulness of quantities. Because 
the aerodynamic forces act in the wind axes, this reference system is used for the translational equations. 
For instance, angle of attack, velocity, and angle of sideslip are either directly measurable or closely related 
to directly measurable quantities, while the body axis velocities (u,v, and w in the a, y, and z directions, 
respectively) are not. The body axis rotational rates are measured by sensors fixed in the body axes; wind 
axis rates can be derived only from these quantities through axis transformations. 


The first reference system to be described is the topodetic reference system, also called the earth-fixed 
reference frame (Etkin, 1972), the earth axes (Thelander, 1965), and the Eulerian axes (Northrop Aircraft, 
1952). The topodetic reference frame is considered fixed in space (and hence, inertial) with the orientation 
of the axes as shown in figure 1; the x axis is directed north, the y axis east, and the z axis down. The 
vehicle position (x and y) and altitude (h) are measured from the origin of this reference system. 


The vehicle-carried vertical axis system (fig. 2; Etkin, 1972) has its origin at the center of gravity of the 
vehicle. The z, axis is directed north, the y, axis east, and the z, axis down. This axis system is obtained 
by a translation of the topodetic axis system to the vehicle center of gravity. The attitude of the aircraft 
(heading, pitch, and bank angles ~, @, and 4, respectively) is described in terms of the orientation of the 
aircraft body axes with respect to the vehicle-carried vertical axes. 


The origin of the body axis system (fig. 3) is the vehicle center of gravity. The z axis is directed toward 
the nose of the aircraft, the y axis toward the right wing, and the z axis toward the bottom of the aircraft. 
The specific orientation of the actual body axes relative to the vehicle body is somewhat arbitrary. For 
symmetrical aircraft, the x and z axes are in the plane of symmetry; for asymmetrical aircraft, these axes 
are located in a plane approximating what would be the plane of symmetry. The positive direction for the 
body axis rates (roll, pitch, and yaw rates, p, g, and r, respectively), the body axis velocities (u,v, and w), 
and the body axis moments (L, M, and N about the z,y, and z axes, respectively) are shown in figure 3. 
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Figure 1. Topodetic axis system. Figure 2. Relationship between topode- 


tic and vehicle-carried vertical axis sys- 
tems. 


N 
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Figure 3, Body axis system. 


The relationship between the vehicle-carried vertical and body axes is shown in figure 4. The Euler 
angles (w, 0, and ) define the orientation of the body axes with respect to the vehicle-carried vertical 
axes. The rotations required to transform the vehicle-carried vertical axes to the body axes are shown in 
figure 5. The heading angle ~ is a rotation about the z vehicle-carried vertical axis into a new axis system 
(designated (21, y1, 21) in fig. 5); the pitch attitude @ is a rotation about the y; axis into the (22, Y2, 22) 
axes system; the roll attitude ¢ is a rotation about the y2 axis into the body axes. 


Xo Xb 
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Figure 4. Relationship between vehicle-carried vertical and body azis systems. 


These rotations are described by 


cosy —sinw 0 


Ly = sinwy cosw 0 (1-5) 
0 0 1 
cos@ 0 sin@ 

Le = 0 1 0 (1-6) 
—sin@ 0 cosé 
1 0 0 

lg = 0 cos? —sind (1-7) 


0 sind cosd 


and the total rotation is described by 


cos 6 cos p cos @ sin 3 —sin 6 
sin@sin@cosw sindsin@siny sin dcosé 
Lay = Ly Lgl, = —cos¢sin wy +cos@cos wv (1-8) 


cosdésinOcosw cos@sin@sinw sin dcosé 
+ sin @sin w — sin dcos w 


or 6) 
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(a) Rotation through > about zy azis. (b) Rotation through @ about y; azis. 


Vy Yo 
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(c) Rotation through ¢ about 2,4 azis. 


Figure 5. Rotation of azes through Euler angles. 


Because Lpy is a unitary matrix, the transformation from the body axes to the vehicle-carried vertical 
axes is bee 

The relationships between the body, wind, and stability axes are shown in figure 6. All three axis 
systems have their origin at the center of gravity of the aircraft. The z axis in the wind reference system 
(tw) is aligned with the velocity vector of the aircraft. The angle of sideslip @ and angle of attack a define 
the orientation of the wind axes with respect to the body axes. (The stability axes are shown in figure 6 
also. This reference system is displaced from the wind axis system by a rotation § and from the body axis 
system by a rotation —a.) 
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Figure 6. Relationship of body, stability, and wind azes. 


Also shown in figure 6 are the components of the velocity vector V in the body axes (u, v, and w) 
and the definition of positive rotations for a and 8. It should be noted that @ is a positive rotation in a 
left-handed coordinate system, whereas the positive sense of all other rotations used in aircraft analysis are 


positive in a right-handed coordinate system. 


The definitions of the body axis velocities (fig. 6) are 


u = Vcosacos# (1-9) 
v=Vsing (1-10) 
w = VsinacosZ (1-11) 


The total velocity V, angle of attack a, and angle of sideslip @ can be expressed in terms of these body axis 
velocities as 


V = |V| = (uv? +0? + wy)? (1-12) 
ww 

a = tan? os (1-13) 
. 1 VU 

6 = sin! a (1-14) 
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1.2 Nonlinear State Equations 


Tor the aircraft problem, the state vector x is 12x 1 vector composed of four 3 x 1 subvectors representing the 
vehicle rotational velocity, the vehicle translational velocity, the vehicle attitude, and the vehicle location: 


x = [xp x2 x3 x7)" (1-15) 
where 
x1 = [pq 7]? (1-16) 
x2 = [V a pg)" (1-17) 
x3 = [¢ 6 x" (1-18) 
x4 = [hz yl" (1-19) 


with x1, X29, X3, and x4 being the rotational velocity, translational velocity, attitude, and position subvectors, 
respectively. The vehicle rotational and translational velocity are defined within the aircraft-fixed axis 
systems. In the formulation of the state used in this report, the vehicle rotations are body axis rates, whereas 
the vehicle velocity terms are stability axis parameters. The vehicle attitude and location parameters are 
earth relative. 


The vector function f, relating the state vector its time derivative, and the control vector to the time 
derivative of the state vector with respect to time, is a 12-dimensional vector function composed of four 
3-dimensional vector subfunctions: 


f[x(t), 3(t), u(t)] = [fF ff fr)" (1-20) 


where f;, f2, fg, and f, are the vector functions that relate the x(t), x(¢), and u(t) vectors to the rotational 
acceleration, translational acceleration, attitude rate, and earth-relative velocity subvectors of x(t). In the 
following sections, each of these subfunctions will be developed separately. The details of the derivation 
of these subfunctions can be found in any of the standard references on aircraft dynamics (Etkin, 1972; 
McRuer and others, 1973; Thelander, 1965). 


1.2.1 Rotational acceleration—The subfunction f, of f from which the rotational acceleration 
terms in the x vector are derived is based on the moment equation 


d 
M= ae (1-21) 
where M is the total moment on the vehicle and H is the total angular momentum of the vehicle. This 
expression can be expanded to 
6 
M= rach +2 x (72) (1-22) 
where 6/6t is the time derivative operator in a moving reference frame (such as the vehicle body axis system) 
and the substitution 
H =12 (1-23) 


has been used to replace the total angular momentum term with the product of the inertia tensor J and 
the rotational velocity vector 0. (The inertia tensor is assumed to be constant with time.) The definition 
of the terms in equation (1-22) follow: 


LL L +Ly7 
M=|/0M/}=|M+ My (1-24) 
“iN N + Ne 
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with L, M, and N being the aerodynamic total moments about the x,y, and z body axes, respectively, and 
Ly, My, and Ny the sums of all power-plant-induced moments; 


Iz =Tyy —Izz 
T= | -Iny Ty —Iyz (1-25) 
—Izz —tyz I, 


where J, Ty, and I, are the moments of inertia about the z, y, and z body axes, respectively, and Ipy, Inz, 
and Jy, are the products of inertia in the x-y, x—-z, and y-z body axis planes, respectively; and 


Q=x,=([p q r]t (1-26) 


where p, g, and r are the rotational rates about the x, y, and z body axes, respectively. Because it is 
assumed that the inertia tensor is a constant with respect to time, equation (1-22) can be rewritten as 


=o =I-'(M—-2 x IQ) (1-27) 


This is the vector subfunction for the rotational acceleration. Designating this subfunction as f,, the 
following definition applies: 


fy[x(t),(¢), u(¢)] = 71M — 9 x (19)] (1-28) 
where 
6 : 
mi” = f,[x(t), X(t), u(t)] (1-29) 
2-4 at (1-80) 


Since the inverse of the inertia tensor J~! is given by 


is , [22 b 
Paarl EE ait 
where 
det I = I, l,l, — In Tj, ~ 1,12, — IyI2, — WLyzTeelry (1-32) 
h = II, - H2, (1-33) 
Ty = Tgyle + yz loz (1-34) 
T= Teylys + Lyle (1-35) 
Ig = Ipl, - 12, (1-36) 
Is = Iglyz + Teylez (1-37) 
Ig = Inly-T, (1-38) 
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the expression for the rotational accelerations can be expanded as a set of scalar equations: 


1 
p= qep feel + UM, +UNI3 - p?(Inrd2 — Ieyls) + pqUe2lt — Iy2T2 — Dal) 


a pr(loy li + DyIg = IyzI3) + P (Lyel a IzyI3) = qr (Deli ez IeyIo + I,213) 
— 7 (lyzly — Tez) (1-39) 


: 1 
q= qa pele + UMI,+ UNIs — P*(Ieels = IgyI5) + pqUrele as Iyzt4 = D_Is) 


— pr(Ipyla + Dyl4 -— IyzTs) + @?(Iyel2 — Inyls) — gr( Delo - Pigdat Teple) 
a r(Tyzl2 pe I,214)] (1-40) 
1 
eS det pels +UMIs + UNIg — p*(IneIs — Inyls) + pq(Ieels — IyeTs - DzIs) 
— prIryls + Dyls — IyzTs) + ?(Iy2ts — Inyle) — qr(Dels - Toyls + Teele) 


— 17(LyzTs — Inels)] (1-41) 
where 
Dr = I, - I, (1-42) 
Dy = Ip - I, (1-43) 
D.=1,-I- (1-44) 


Equation (1-3) defines the generalized nonlinear state equations as 
TR(t) = £[x(4), X(2), u(@)] 


This equation, although more complicated than the nonlinear equations defined by equation (1-1), allows for 
a more tractable formulation of the state equation by using the matrix T to provide a means of addressing 
the rotational accelerations in a decoupled axis system. 


The derivation of the rotational acceleration terms is based on the moment equation ( 1-22): 
M = (IM) +9 x 10 


Rearranging terms and assuming that the inertia tensor is constant with respect to time, the equation can 
be written as 5 


1-0 =M~2x 19 (1-45) 


The rows of this vector equation are now scaled using the following scaling matrix: 


l/r 0 0 
=|] 0 I 0 (1-46) 
GO Or “Ge 


This matrix, when premultiplying equation (1-27), merely divides the first row by the roll inertia J,, the 
second row by the pitch inertia J,, and the third row by the yaw inertia J,. Using the definition 


J=J'I (1-47) 
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the resulting equation is 


6 


Joa = J’M — J'(Q x I) (1-48) 
and J can be written as 
1.0 —TIpyf/Ie  — Iez/Ix 
J=| —Izy,/I, 1.0 ~Tyzfly (1-49) 


— I,2/Iz — ye [Tz 1.0 


Equation (1-48) can be expanded and expressed as 


p 1.0 oar Tey / Tp — Ine] Tn Dp 
q | = Li, 1.0 —Tyefl, |) | ¢ 
rt _ ee - ye/de 1.0 t 


UL/Iz = TPL zy / Ie + PLezf Iz + rqly| Ix + (¢? = F lysis _ qrd,/Iy 
= UM/Iy — rple/Ty + 7qTey/Ty — pel ye/Iy + (r? — p*)Iez/Ty + prle/Ty (1-50) 
IN/I, + QPle [Iz — QtIz2/I2 + prlyz/Tz + (p? = @ Iey/ I, on pay / Tz 
where ’, q’, and 7 are the decoupled rotational accelerations of the vehicle. 


Using the definition of J in equation (1-49), the matrix transformation T can be defined as 


I 
J 103x3 | 
Soeeinee O6x6 
fied acs ho ee (1-51) 
t 
exe lexe 
i 
\ 


which would be an identity matrix except for the presence of the inertia terms in the upper left-hand corner. 
Thus, the vector subfunctions for the generalized state equation defining vehicle translational acceleration, 
vehicle attitude rates, and earth-relative velocities are the same as those defined for the standard nonlinear 
state equations in sections 1.2.2, 1.2.3, and 1.2.4, respectively. 


1.2.2 Translational acceleration.—Derivation of the translational acceleration vector subfunction 
f, is based on the force equation 


F= £ (mv) (1-52) 
where F is the total force acting on the vehicle and m is the vehicle mass. This expression can be expanded to 
Fem (-v +Qx v) (1-53) 

with the assumption of constant mass with respect to time and the following definitions of F and V: 
F=[5X SY yz}? (1-54) 


where UX, XY, and LZ are the sums of the aerodynamic, thrust, and gravitational forces in the z,y, and 
2 body axes, respectively, and 
Value wl? (1-55) 
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Rearranging the terms of equation (1-52) gives an expression for the translational acceleration: 


Oa Se Oey (1-56) 
ét m 


This equation expresses body axis accelerations in terms of body axis forces, angular rates, and velocities. 
However, the desired form of this relation requires the translational accelerations in the wind axis system; 
that is, in terms of the magnitude of the total vehicle velocity V, angle of attack a, and angle of sideslip 6, 
which are expressed by equations (1-9) to (1-11) 

u = Vcosacosf 

v= Vsing 

w= Vsinacos?Z 


and equations (1-12) to (1-14) 


< 
I! 


[Vv | = (u? + v2 4: w?)i/? 


w 
a = tan? (=) 
U 


8 = sin“! (+) 


The wind axis translational acceleration terms (derived in app. B) are summarized as: 
[V & B]* = fh[x(t), x(t), u(t)] (1-57) 
where 


. 1 
Vo= ml DcosB+Y¥ sing + Xpcosacos# + Yrsin8 + Zrsinacos 8 


— mg(cos a cos § sin # — sin @ sin ¢ cos 6 — sin a cos J cos ¢ cos 6)| (1-58) 
: 1 : Wnts 
a= Varcsce: -~L+ Zy cosa — Xzsin a+ mg(cosacos ¢cos 6 + sin asin 6)] 
+q-tan 6 (pcosa+rsina) (1-59) 


: 1 
B= ay lD sin 8 + ¥ cos 8 — Xqzcosasin§+Ypcos@ — Zrsinasin B 
+ mg(cos asin sin @ + cos B sin ¢cos 6 — sinasin B cos ¢cos @)| + psina —rcosa (1-60) 


with D being total aerodynamic drag; Y total aerodynamic sideforce; and Xt, Yp, and Zr total thrust 
force along the x, y, and z body axes, respectively. 


1.2.3 Attitude rates.—The matrix R that transforms angular velocities in the earth-fixed axis system 
into body axis angular velocities is defined by 


1 0 —sin@ 
R=1]0 cos¢d = sin¢dcosé (1-61) 
Q -sing cosdécosé 


where £ is derived by Maine and Iliff (1986) from the total angular velocity of the aircraft expressed in 
terms of the derivatives with respect to time of the Euler angles (¢, 6, #): 


7) db 1 0 0 0 1 O 0 cos6 0 —siné | | 0 
q| = {0/+ {0 cos¢d sing O|/+ |0 cosd sing 0 1 =O 0 
r 0 0 —sing cosd| | 0 0 -sing cos¢ sin@ 0 cos@ cy) 
1 O — sin @ é 
= |0 cos¢@ sin¢cosé | | @ (1-62) 
0 —sin @ cos gcos@ o 
This transformation from earth-fixed to body axes can be expressed by the equation 
o=n(4e) (4 
~ dt 
where E is an attitude vector whose components are the Euler angles: 
=[¢oy]" (1-64) 


Premultiplying both sides of equation (1-63) by R71 and rearranging terms yields the equation for the 
attitude rates, 


{B= RO (1-65) 
which can be expanded into the scalar equations 
¢ = ptqsingtand +rcosdtand | (1-66) 
6 = qcosf—rsing (1-67) 
b = qsindsecé + rcos¢secé (1-68) 


1.2.4 Earth-relative velocity——The matrix Lpy that transforms earth axis system vectors into the 
body axis system is defined by equation (1-8) as 


cosy —siny 0 cos@ QQ sin8@ 1 0 0 
Ipy = singw cosy 0 0 1 O 0 cosd —sing 
0 0 1 —sind 0 cosé 0 sing cos¢ 
cos 6 cosy cos @sin wy —sin#@ 


sin gsin@ cos —cosdsiny sin dsin 6 sin y + cosdcos yp sin dcos 6 
cosdsin@cosi+sindsiny cosdsin Osin yy —sing@cosy cosdcosd 


The specific relationship between earth-relative velocities and body axis velocities is expressed by 
d 
V= Ivy —R (1-69) 
where R is the earth axis system vector defining the location of the vehicle: 
=[z y 2]" (1-70) 
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with z= —A. 
The equation for the earth-relative velocity can be formulated as 


d -1 
ae = L5yV (1-71) 
in which these velocities are expressed in terms of body axis velocities. Using equation (1-72) and the 
definitions of the body axis velocities in equations (1-12) to (1-14) allows the earth-relative velocities to be 
expressed in terms of V, a, and : 


h= V(cos a cos @ sin 9 — sin @ sin ¢ cos @ — sin a cos f cos ¢ cos @) (1-72) 


z = V[cosecos f cos 6 cos p + sin B(sin ¢sin 6 cos p — cos ¢sin w) 
+ sin a cos (cos sin @ cosy + sin Pd sin ~)] (1-73) 


y = V[cosacos f cos 6 sin w + sin G(cos pcos » + sin dsin 6 sin p) 
+ sin a cos B(cos ¢ sin # sin p — sin ¢ cos Y)] (1-74) 


1.3. Nonlinear Observation Equations 


No standard set of observation variables exists for the aircraft analysis and control design problem. However, 
for any guidance and control problem, the main observation variables generally will be a subset of the state 
variables. Other common observation variables are the vehicle body axis translational accelerations and 
air data parameters. Thus, the dimension of g[x(t), x(¢), u(t)] is not fixed and varies from application to 
application. The set of observation variables described in this section was selected to address a wide range 
of problems. The basic composition of the observation vector y as used in this report is given by 


y= [x? xt ut y|F (1-75) 


where x and x are the state vector and time derivative of the state vector described previously, u is the 
control vector, and y’ is defined by 


oy sly ys ys yf ys ye y? ysl” (1-76) 
where 

yi [Geka Gee gaye ey ag) yma oa) (1-77) 
y, = [a M Re Re! G d Qc/Pa Pa BT Ti)" (1-78) 
y3 = [y fpa h]* (1-79) 
ya = [Es P,]* (1-80) 
ys=[L DN A]T (1-81) 
ye=[uvwiad af? (1-82) 
yz = les Be ba byl? (1-83) 
yg =[T pe a ]* (1-84) 
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with the elements of y{ being terms related to the vehicle body axis acceleration, the elements of y}, being 
air data terms, the elements of y being flightpath-related terms, the elements of y4, being terms related to 
vehicle energy, ys being a vehicle force vector, the elements of y, being body axis translational rates and the 
time derivatives of those terms, y> being a vector of variables representing measurements from instruments 
not located at the vehicle center of gravity, and the elements of y$ being a collection of miscellaneous terms. 
Obviously, this grouping of terms is somewhat arbitrary and is done primarily to ease the definition of these 
terms in the following sections of this report. This grouping of observation variables parallels that used by 
Duke and others (1987). 


The vector function g relating the state vector, the time derivative of the state vector, and the control 
vector to the observation vector is an é-dimensional function composed of four subfunctions: 


glx(t), (2), u(t)] = [x? x7 ut gi] (1-85) 


where x, X, and u are identity functions on the state vector, time derivative of the state vector, and control 
vector, respectively, and g’ is composed of vector subfunctions defining the y’ vector. 


The state vector, time derivative of state vector, and control vector components of the observation 
vector are not discussed in detail in this section of the report. The equations for the elements of the time 
derivative of the state vector were developed in section 1.1. The observation equations for the state and 
control variables are simply identities." The equations for the remaining observation variables are obtained 
from a variety of sources. In addition to the previously cited sources, Clancy (1975), Dommasch and 
others (1967), Gainer and Hoffman (1972), and Gracey (1980) provide the background and derivation of 
the observation equations used in this report. 


1.3.1 Accelerations.—The vehicle body axis accelerations and accelerometer outputs constitute the 
set of observation variables that, after the state variables themselves, are most important in the aircraft 
control analysis and design problem. These accelerations and accelerometer outputs are measured in units 
of g and are derived directly from the body axis forces defined in section 1.2.2. The body axis acceleration 
vector a can be expressed as 


d 6 
ae Hav TV ; (1-86) 
It is important to note here that the u, #, and w body axis velocity rates, derived in appendix B and 
defined by equation (B-1), are not the body axis accelerations. The body axis accelerations contain not 
only the body axis velocity rates but also the rotational velocity and translational velocity cross-product 


terms. Thus, expanding equation (1-86) yields 


az k u+qwu—rv 
a=] ay. | =| o+ru—pw (1-87) 
azk w+ pv — qu 


where dzk, @y,k, and a;, are the kinematic accelerations in the vehicle body z, y, and z axes, respectively. 
Using 


a (1/m)(X7 + Xa + X_) + rv —- qu 
6] =| (/m)\(¥r+¥.+¥,) + pw-ru (1-88) 
w (1/m)(Zy + Z,+ Zg) + qu — po 
(stated as eq. (B-1) in app. B), equation (1-87) can be rewritten as 
ark (1/m)(Xr + Xa + Xz) 
yk | =| (1/m)(¥r+Y¥, + ¥,) (1-89) 
azk (1/m)(Zr + Za + Zg) 
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where X,, Ya, and Z, are total aerodynamic forces and X,, Y,, and Z, are total gravitational forces along 
the z, y, and z body axes, respectively. This can be expanded in terms of the gravitational and aerodynamic 
forces to give (in units of g) 


azk 1 Xz — Deosa+ Lsina — gmsin# 
tye: = Yr+Y +gmsin ¢cos 6 (1-90) 
Gzk Jo Zz — Dsina — Lcosa + gmcos¢cosé 


where go is the acceleration due to gravity at sea level. 


The outputs of body axis accelerometers at the vehicle center of gravity are simply the body axis 
accelerations due to the thrust and aerodynamic forces. The accelerometer output equations can be written 
directly from equation (1-90) as 


ay l X7—- Dcosat+Lsina 
ay | =—— Yr+Y (1-91) 
a, gom Zy — Dsina — Lcosa 


where az, dy, and a, are the outputs of accelerometers at the vehicle center of gravity and aligned with the 
vehicle body z, y, and z axes, respectively. Because the normal acceleration a, is defined by 


an = ~Gz (1-92) 

an expression for this variable can be extracted from equation (1-91): 
ay = (-—2Z7 + Dsina+ Lcosa)/ggm (1-93) 
The equations defining the output of accelerometers aligned with the vehicle body axes but displaced 


from the vehicle center of gravity are derived by Gainer and Hoffman (1972) using the definition of inertial 
acceleration given in equation (1-86) 


siya oy 


~ bt 
and the definition of inertial velocity 
6 
Va grtQxr (1-94) 
The results from Gainer and Hoffman (1972) are reproduced here without rederivation: 
as, de —[(q? + 9)te — (pq — F)y¥x — (pr + 4)22]/90 
dy,i | = | dy + [(pg t+ Fey — (p? +1? )yy — (gr — p)2y]/gG0 (1-95) 
az,i az + [(pr — 4)ez + (ar + Bus — (¢? + P*)22)/90 


where az,j, 4y,;, and @,,; are outputs at accelerometers aligned with the x, y, and z body axes but not located 
at the vehicle center of gravity; the subscripts z, y, and z refer to the z, y, and z body axes, respectively; 
and the symbols z, y, and z refer to the x, y, and z body axis locations of the sensors relative to the vehicle 
center of gravity. Because the normal acceleration is the negative of the z body axis accelerometer, the 
output of a normal accelerometer not at the vehicle center of gravity but aligned with the z body axis, ani, 
is given by 

On, = On — [(pr — 4)az + (qr + Pye — (@ + P*) 221/90 (1-96) 


The final quantity included in the general category of accelerations is load factor n. This quantity is 
defined without inclusion of the z body axis force component as 


= (1-97) 
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1.3.2 Air data parameters.—The air data parameters having the greatest application to aircraft 
dynamics and control problems are the sensed parameters and the reference and scaling parameters. Chosen 
for inclusion as the sensed parameters are impact pressure qc, static or free-stream pressure pg, total pressure 
pt, ambient or free-stream temperature T, and total temperature T;. The selected reference and scaling 
parameters are Mach number M, dynamic pressure g, speed of sound a, Reynolds number Re, Reynolds 
number per unit length Re’, and the Mach meter calibration ratio g./p,. The derivation of these quantities 
is treated extensively by Gracey (1980). 


The nonlinear equations defining these quantities are 


1/2 
Po 
= 142-7] 1-98 
a | Polo ie) 
V 
M=— (1-99) 
Rees (1-100) 
av 
Re (1-101) 
wv 
Ps os, 
q= eV (1-102) 
_ J (G0 + 0.2M?)35 — 1.0]p, (M < 1.0) (1-103) 
fe = | (1.2M?(5.76M?/(5.6M? — 0.8)]?5—1.0}p, (M > 1.0) 
a _ } (1.04 0.2M?)*5 — 1.0 (M < 1.0) (1-104) 
Pa | 1.2M*[5.76M?/(5.6M? — 0.8)]?5—1.0 (M > 1.0) ; 
T, = T(1.0+ 0.2M”) (1-105) 


where p is the density of the air, 4 is the coefficient of viscosity, and the subscript 0 refers to sea level, 
standard day conditions. Free-stream pressure, free-stream temperature, and the coefficient of viscosity are 
properties of the atmosphere and are assumed to be functions of altitude alone. 


1.3.3. Flightpath-related parameters.—Included in the observation variables are what might best 
be termed flightpath-related parameters for lack of better nomenclature. These terms include flightpath 
angle y, flightpath acceleration fpa, and vertical acceleration h. The variables are defined by the following 


equations: 
h 
= sin7l : 
= sin (F) (1-106) 


¥ 
Vv 

fpa = — (1-107) 
g 

h= Oz 8in 9 — ay, sin dcos 8 — az, cosgcosd (1-108) 
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1.3.4 Energy-related parameters.——Two energy-related parameters are included with the observa- 
tion variables considered in this report: specific energy E,, and specific power P,, defined as 


y 
= ae 1A 
E, = h+ 7 (1-109) 

dE, =: VV 
mee pais -110 
B= Shs ; (1-110) 


1.3.5 Force parameters.—The set of observation variables being considered also includes four force 
parameters. These quantities are total aerodynamic lift L, total aerodynamic drag D, total aerodynamic 
normal force N, and total aerodynamic axial force A, defined as 


L=GSCy (1-111) 
D = GSCp (1-112) 
N = Leosa+ Dsina (1-113) 
A= -Lsina+ Dcosa (1-114) 


where S is the surface area of the wing, C\, coefficient of lift, and Cp coefficient of drag, 


1.3.6 Body axis rates and accelerations.—Because they are of interest in the control analysis and 
design problem, six body axis rates and accelerations are included as observation variables. These include 
the x body axis rate u, the y body axis rate v, and the z body axis rate w.. Also included are the time 
derivatives of these quantities, a, 0, and w, respectively. 


The definitions of the body axis rates are given in equations (1-9) to (1-11) as 
u = Vcosacos@ 
v=Vsinf 


w= Vsinacos? 


The time derivatives of these terms can be defined using equation (B-1) and equations (B-8), (B-9), (B-10), 
and (1-56) as 


._ Atr-—gmsind — Deosa+ Lsina 


iS Te ry Sie = GV sine coe e (1-115) 
i Y 
= Ar taming cos 0+ Y 5 pV sina.cos f — rV cos cos 8 (1-116) 
We D sie 
ji As OOS ADE OS DOES eh aise de a Cee (1-117) 


m 


1.3.7 Instruments displaced from the vehicle center of gravity.—The need to include measure- 
ments from instruments displaced from the vehicle center of gravity arises from the fact that not all aircraft 
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instrumentation is located at the vehicle center of gravity. The most important of these quantities are un- 
doubtedly the accelerometer outputs treated in section 1.3.1. In this section four additional parameters are 
presented: angle of attack (a,;), angle of sideslip (G,;), altitude (h;), and altitude rate (h.;) measurements 
from instruments displaced from center of gravity by some z, y, and z body axis distances. The subscripts 
a, B, h, and A refer to the displacements of the angle-of-attack, angle-of-sideslip, altitude, and altitude rate 
instruments from the vehicle center of gravity. The equations used to compute these quantities are 


pene Hee (1-118) 
Vv 

pce peo (1-119) 
V 

h; = h+2,sin@— yp, sin dcos@ — z, cos dcosd (1-120) 


hi=ht A(x; cos6 + y; sin dsin # + 2; cos dsin 6) — HY; cos pcos @ — z; sin ¢ cos @) (1-121) 


1.3.8 Miscellaneous observation parameters.—The final set of observation parameters considered 
in this report is a miscellaneous collection of parameters of interest in analysis and design problems. These 
parameters are total angular momentum T, stability axis roll rate pg, stability axis pitch rate q,, and stability 
axis yaw rate rs. The equations used to define these quantities are , 


T = 5le? — QIeypq ~ Weepr + Tyg? — 2lyzqr + Ler?) (1-122) 
Ps = pcosa+rsina (1-123) 
Is = 4 (1-124) 
e = —psina+rcosa (1-125) 


2 LINEAR SYSTEM EQUATIONS 


The standard state equation for a linear differential system has the form 
x(t) = A’x(t) + B’u(t) (2-1) 


where, for a time-invariant system, A’ is a constant n X n matrix and B’ is a constant n x k matrix. The 
standard output equation has the form 


y(t) = H’x(t) + F’u(t) (2-2) 


where H’ is a constant £ x n matrix and F” is a constant £ x k matrix. The generalized linear system 
equations used with an extended formulation compatible with the generalized nonlinear equations (1-3) and 
(1-4) can be characterized by 


Cx(t) = Ax(t) + Bu(t) (2-3) 

y(t) = Hx(t)+ Gk(t) + Fu(t) (2-4) 

where C and A are constant n xX n matrices, B is a constant n x k matrix, 1 and G are constant 2x n 
matrices, and F is a constant £x k matrix. The nonlinear system equations developed in section 1 (eqs. (1-1) 


to (1-4)) can be linearized about a trajectory, and a linear model can be formulated that is similar to either 
the standard or the generalized linear system equations. 
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2.1 Linearization of the State Equation 


If uo(¢) is given input to a system described by the state differential equation (1-3), and if xo(¢) is a known 
solution of the state differential equation, then approximations to the neighboring solutions can be found 
for small deviations in the initial state and in the input by using a linear state differential equation. The 
nonlinear state differential equation (1-3) can be linearized about a general trajectory, as by Kwakernaak 
and Sivan (1972) and Dieudonne (1978), so that xo(t) satisfies 


TXo(t) = £[Xo(t), xo(t), uo(Z)] 


Assuming that the system is operated at close to nominal conditions with u(t), x(t), and x(t) deviating 
only slightly from uo(t), xo(t), and Xo(¢), the following expressions can be written: 


u(t) = u(t) + du(t) (2-5) 
x(t) = xo(t) + 6x(t) (2-6) 
x(t) = Xo(t) + 6x(t) (2-7) 


where du(t), 6x(t), and 6x(¢) are small perturbations to the control, state, and time derivative of the state 
vectors, respectively. 


Substituting equations (2-5) to (2-7) into the nonlinear state differential equation (1-3), expanding in a 
Taylor series about xXo(t), Xo(t), uo(t), and assuming T constant with respect to x(t) yields 
: : i of of .. Of 
t = pad pas or ; 
T[Xo(t) + 6x(t)] = f[xo(t), Xo(t), u(t)] + Ix 6x+ Dx 6x+ au 6u + h(t) (2-8) 
where 0f/0x, Of/Ox, and Of/Ou are defined in equations (2-9) to (2-11) and h(t) represents the sum of 
the higher order terms in the Taylor series, assumed to be small with respect to the perturbations. The 
matrices used in the Taylor series expansion are defined by the following relationships: 


of Of 
ax = axl. (2-9) 
(X0,X0,Uo) 
of af 
Be = El ay 
(Xo .Xo0 Uo) 
of Of 
Ba = au (2-11) 
(X0,X0,Uo) 
the (i, 7)th elements of which are defined as 
of Of; 
ae =e (2-12) 
(5) — Of 2-13 
Ox) 5; OX; ee) 
of Of; 
all ae (2-14) 


respectively, where f; is the ith simultaneous equation of the nonlinear state differential function in equa- 
tion (1-3), x; the jth element of the state vector, x; the jth element of the time derivative of the state 


23 


vector, u; the jth element of the control vector, and all derivatives are evaluated at the nominal condition 
(xo(t), Xo(t), uo(t)). 

Subtracting equation (1-3) from (2-8), rearranging terms and neglecting the higher order terms yields a 
linearized state equation, 

Ir 7 sr IG oe a eee (2-15) 
Ox Ox Ou 

where the arguments of the matrix functions have been dropped to simplify the notation and where it is 
understood that the matrices are to be evaluated along the nominal trajectory. 


Letting 
of 
C=T-— (2-16) 
of 
A= a (2-17) 
of : 
B= Ba (2-18) 
equation (2-15) can be written as 
C 6x(t) = A 6x(t) + B bu(t) (2-19) 


which is precisely the formulation of the generalized state equation desired. 


Premultiplying both sides of equation (2-19) by C~} results in the standard form of the linearized state 
differential equation, 


6x(t) = C-1A 6x(t) + C71 B 6u(t) (2-20) 

Letting 
A'= CA (2-21) 
B'=C™'B (2-22) 


equation (2-20) can be written in the more usual notation 


6X(t) = A’ 6x(t) + B’ éu(t) (2-23) 


2.2 Linearization of the Observation Equation 


The technique used in section 2.1 to linearize the state equations can be applied to the nonlinear observation 
equation (1-4), 
y(t) = g[x(t), x(t), u(t)] 


Performing a Taylor series expansion about the nominal trajectory (xo(t), xo(t), uo(t)) yields 


; 0 og .. 9 
yo(t) + 6y(t) = g[xo(t), Xo(é), wo(t)] + 5 ox + SE x + 5% du + h(t) (2-24) 
where 
og _ 0 


x . 
(Xo,X0,Uo0) 
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oe . 9B (2-26) 
Ox Ox ie teas) 
8 - % (2-27) 
= "I (x0 ko 40) 
the (7, 7)th elements of which are defined by 
o Og: 
2 = 2-28 
ie) ig Ox j ( ) 
dg agi 
a ee 2-29 
coe dx; (29) 
dg ogi 
p=} = 2-30 
Car Ou; ( ) 


respectively, where g; is the ith simultaneous equation of the nonlinear observation equation (1-4). Again, 
all derivatives are evaluated at the nominal condition (x9(t), Xo(t), uo(t)). 


Subtracting equation (1-4) from equation (2-24), rearranging terms, and neglecting higher order terms 
results in a linear observation equation, 


_ Og Og .. Og 
where the arguments of the matrix functions have been dropped to simplify notation. Letting 
og 
H= Dx (2-32) 
dg 
G= ax (2-33) 
dg 
Fo au (2-34) 
equation (2-31) can be rewritten as 
dy(t) = H éx(t) + G éx(t) + F éu(t) (2-35) 


which is the generalized linear observation equation desired. 


The standard form of the observation equation can be derived by substituting for 6x from equation (2-23) 
into equation (2-33). This substitution results in 


6y(t) = H 6x(t) + G[A' 6x(t) + B’ éu(t)] + F 6u(t) (2-36) 
which can be written as 
éy(t) = [H + GA’) 6x(t) + [F + GB’] éu(t) (2-37) 
By letting 
H'’=H+GA' (2-38) 
Fi = F+GB' (2-39) 


equation (2-37) becomes 
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6y(t) = H’' éx(t) + F’ 6u(t) (2-40) 


2.3 Definition of Matrices in Linearized System Equations 


The results of sections 2.1 and 2.2 can be used to define the matrices in the linearized system equations 
in terms of partial derivatives of the nonlinear state and observation functions taken with respect to the 


state, time derivative of state, and control vectors. All derivatives are understood to be evaluated along the 
nominal trajectory. 


Using the nonlinear state equation (1-3), 
Tx(t) = f[x(t), x(t), u(2)] 
the terms in the generalized form of the linearized state equation (2-19), 
C 6x(t) = A &x(t) + B éu(t) 


can be defined as 


c=r-ot | (2-41) 
f 

A= x (2-42) 
f 

B= es (2-43) 


The terms in the standard form of the linearized state equation (2-20), 
6x(t) = A’ 6x(t) + B’ bu(t) 


can be defined as 


; Of)—! of 
es [r- Rl Bx 2%) 
of)" of 
p= [0-5] oe (2-45) 


In a similar manner, the nonlinear observation equation (1-4), 
y(t) = glx(t), x(t), u(t)] 
can be used to define the terms of the generalized linearized observation equation (2-35), 


by(t) = H 6x(t) + G 6X(t) + F bu(t) 


as 
H=~> (2-46) 
G=s (2-47) 
F== (2-48) 
The terms in the standard form of the linearized observation equation (2-40), 


éy(t) = H’ x(t) + F’ u(t) 


can be defined as 


1_ 98 | og a = es 4 
a= Ox + Ox - Ox ox oe) 
Og Og ai of 
U = e_—_ — — es - 
prs du = Ox @x} au aa) 


2.4 Elements of the Linearized System Matrices 


The elements of the linearized system matrices derived in sections 2.1 and 2.2 are determined by applying the 
linearization method employed with the vector equations in those sections to the individual scalar equations 
constituting the vector equations that define the time derivatives of the state and observation variables. 
Thus, for a matrix, such as the state matrix A defined by equation (2-42), 


_ of 


eo 


the element occupying the ith row and jth column of A, (A)j,;, can be represented as 


Of; 
(A)ig = Ox; (2-51) 


where f; is the scalar function defining the time derivative of the ith state and x; is the jth state. The 
individual terms used in the A, B, C, H, G, and F matrices are defined in appendix D based on the 
generalized derivatives derived in appendix C. 


Using the state vector x defined in (1-7) as 


x=[pqrVaBhoovhay)" 


the elements of the A matrix can be expressed as 


A(p')/Op O(p!)/Oq +++ O(p!)/By 
O(q')/Op O(q')/0q «++ O(q')/Oy 
A=| : (2-52) 
O(X)/Op O(x)/Oq +++ O(x)/Oy 
A(y)/Op O(y)/Aq --- A(y)/Ay 
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Substituting for these partial derivatives using the terms in appendix D gives 


(1/Ix)[(GSb2/2V0)Ce, + OLr/Op (1/Le)[(GSb2/2V0)Ce, + OL-r/y + Ine Po 


—InyPo + Iz290] +2Iy290 + ro(Ly = I,)] 
A= | (1/Iy)[(GSbe/2Vo)Cm, +OMr/Op —(1/Ty)(GS22/2Vo)Cm, + OMr/dq ae (2-53) 


—2Ip2 Po — Lyz40 + ro( I, =; Tz)} +InyTo = LyzPo] 


The elements of the B, C, H, G, and F matrices can be determined in a similar fashion, although some 
care must be taken in determining the elements of the matrices for the observation equation and the C 
matrix. 


To determine the elements of the matrices for the observation equation, one must consider the definition 
of the nonlinear vector function g defining the observation variables (eq. (1-85)), 


glx(t), X(t), u(t)] = [xT 7 ug!) 
and the definitions of the matrices for the generalized linear observation equations (2-46) to (2-48) , 


H= 
G=8 


Fos 


These matrices may be expressed using a partitioning based on the vector subfunctions of g as 


H=|_. (2-54) 


(2-55) 


(2-56) 
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which become 


li2x12 


O12x12 
on (2-57) 
Oxx12 


og 
x 


012x12 


cen eee (2-58) 


Pal Sea (2-59) 


upon evaluating the partial derivatives of the identity functions x, x, and u. 


The C matrix may be viewed as a partitioned matrix as 


] ] 
Cyr! Cie | 

ere ae Dexe 
1 

C Bical eae! (2-60) 

i} 

Ooxe | loxe 
1 
i} 
1 


where, from equation (1-48), 


1.0 = ry/ Le —Iyz/Te 


Ci =J=|- az/ly 1.0 oe az/Ly (2-61) 
—Ine/T, —Iy:/I, 1.0 
and 
—A(p')/OV —A(p')/d& —8(p')/8B 0 —(GSbé/2VoIz)Ce, —(GS0?/2Volz)Ce, 
Cy = -a(q/)/aV -A(q')/06 -8(q/)/08 = {0 —(qSe?/2VoIy)Cme —(GSbE/2Voly)Cm, (2-62) 
—A(r\JOV —A(r) da —A(r)/B 0 —(GSbe/2Volz)Cng —(G50?/2VoIz)Cn, 
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1.0-AV)/AV  ~A(V)/da~— -A(V)/aB 
Coo = -8(&)/OV 1.0— 6(4)/8& — -8(&)/98 
—0(8)/AV —0(8)/O& 1.0 — A(8)/d8 
1.0 (GS@/2Vom)(cos Bo Cp, — sin Bo Cy,) (G5b/2Vom)(cos fo Cp;) 
0 1.04 (¢St/2V3em cos Bo)CL; (gSb/2V¢m cos Bo)CL, (2-63) 
0 (GSe/2VPm)(sin Bo Cp, + cos Bo Cy,) 1.0 — (GSb/2VZm)(sin Bo Cp, + c08 Bo Cy,) 


The inverse of the C matrix, C~!, can be expressed as a partitioned matrix in terms of the matrix subpar- 


titions of the C matrix as 


cee eee Mescstal, (2-64 


The elements of the A’, B’, H', and F’ matrices can be determined using the C-? matrix defined in 
equation (2-64), the A, B, H, G, and F matrices, and the definitions for A’, B', H’, and F" given in 
equations (2-21), (2-22), (2-38), and (2-39). 


3 CONCLUDING REMARKS 


This report derives and defines a set of linearized system matrices for a rigid aircraft of constant mass, flying 
in a stationary atmosphere over a flat, nonrotating earth. Both generalized and standard linear system 
equations are derived from nonlinear six-degree-of-freedom equations of motion and a large collection of 
nonlinear observation (measurement) equations. 


This derivation of a linear model is general and makes no assumptions on either the reference (nominal) 
trajectory about which the model is linearized or the symmetry of the vehicle mass and aerodynamic 
properties. 


Ames Research Center 

Dryden Flight Research Facility 

National Aeronautics and Space Administration 
Edwards, California, January 8, 1987 
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APPENDIX A—AERODYNAMIC FORCES AND MOMENTS 


The aerodynamic forces and moments acting on an aircraft are the result of multiple factors whose signif- 
icance varies with flight condition as well as from vehicle to vehicle. In general, these forces and moments 
are nonlinear functions primarily of Mach number, angle of attack, angle of sideslip, altitude, rotational 
rates, and control-surface deflections. For the purposes of this report, the aerodynamic forces and moments 
are assumed to be functions having the following form: 


F = (a, 8,V,h,p, 9,7, & B, 51,---,6n) (A-1) 


where F is an arbitrary force or moment, ® is an arbitrary function, and the 6; are the n control surface 
deflections. These forces and moments are related to the nondimensional force and moment coefficients by 
the equations for the forces, 


D = @SCp (A-2) 

Y = gSCy (A-3) 

L = GSCy (A-4) 
and the moments, 

L = qSbC, (A-5) 

M = GS@'m (A-6) 

N = qSbC, (A-7) 


where b is reference span and @ is reference aerodynamic chord. 


While the nondimensional aerodynamic force and moment coefficients are themselves nonlinear func- 
tions of the vehicle states, time derivatives of the vehicle states, and the control surface deflections, these 
coefficients are commonly expressed in linear form in terms of partial derivatives of these coefficients with 
respect to the functional variables. These linear equations for the aerodynamic force and moment coeffi- 
cients are derived in the same way as the linearized system equations (section 2); therefore, this derivation 
will not be repeated here. These linear equations are 


Ch = Chip + Cia + Crgh + CLzh t+ Cry V 


+ = C15, 6i + CLyB + CLad + CL, + Cig + C46 (A-8) 
Cp = Coo + Coa + Co,8 + Co, h + CoyV 
+ 3 Cp,.6: + Copb + Cp.4 + Co, * + Co,& + Cp, 8 (A-9) 
i=1 
Cy = Cy, + Cy,a+ Cy,6 + Cy,h + CyyV 
435 Ov, 5: + Cv,B + Crgd + Cv. + Cyd + Cv, 8 (A-10) 
i=l 
Ce = Ce + Cog + Cog + Cok + CoV 
4 3500, f: + Cob + Cod + CoP + Cand + Ce,8 (A-11) 


i=1 
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Cn On PO On oC, ECV 


+ 37 Cig, 5: + Cmpb + Cig G + Cre? + Cmad + Cm, 8 (A-12) 


wl 


Ce Cp POE D POpRL CV 


i 7 
+7 C14, 5 + Cngh + CngG + Cn? + CraG + Cn ,8 (A-13) 


t=1 
where Cg, is the value of the coefficient along the nominal trajectory and the notation Cz, is defined as 


nee 
& = “Oa 


with Cg being an arbitrary force or moment coefficient and a being an arbitrary state, time derivative of state, 
or control-related parameter that for the usual derivatives is nondimensional. However, the derivatives with 
respect to altitude and velocity are not taken with respect to a nondimensional quantity. The definitions 
of these nondimensional stability and control derivatives are given in terms of the coefficient Cz. The 
nondimensional stability derivatives are defined as 


(A-14) 


Cr, = os (A-15) 
Ces = = (A-16) 
C= Tata (A-17) 
Ce, = Fenty (A-18) 
C. = amTD (A-19) 
Ce, = yeaa (A-20) 
Ce, = Tea (A-21) 
The two other stability derivatives are not nondimensional and are defined as 
Ce, = ot (A-22) 
Cz, = oe (A-23) 
The control derivatives are defined as 
Ce, = ai (A-24) 
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The rotational terms in equations (A-8) to (A-13) are nondimensional versions of the corresponding vari- 
able with 


p= 3 (A-25) 
i= (A-26) 
pe (A-27) 
a= ae (A-28) 
p= oe (A-29) 


Because the C¢, terms are included, the force and moment coefficients are total force and moment coefficients. 
The state, time derivative of state, and control parameters on the right-hand side of equations (A-8) to 
(A-13) are differentials. 
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APPENDIX B—-DERIVATION OF THE WIND AXIS 
TRANSLATIONAL PARAMETERS V, a, AND 6 


The derivation of the wind axis translational acceleration parameters is based primarily on the definitions 
in equations (1-9) to (1-14), the body axis translational acceleration equations (1-56), and the expression 
of the force terms defined in equation (1-53). In the following sections, each of the wind axis transla- 
tional acceleration terms is derived separately after stating some preliminary definitions applicable to all 
calculations. 


B.1. Preliminary Definitions 


Equation (1-56), 
ft) 
—=V= if ~-Qx%V 
ét m 

can be expanded, using equations (1-54), (1-55), and (1-26), to 


ts (1/m)(X7 + Xa+X,) + rv — qu 
o{ =| (/m)(Yr+ Yat ¥,)+ pw-ru (B-1) 
w (1/m)(Zr + Za + Zg) + qu — pv 


The body axis aerodynamic forces can be rewritten in terms of the stability axis forces lift L, drag D, and 
sideforce Y: 


X, = —Dcosa+Lsina (B-2) 
Y,=¥ (B-3) 
Z, = —Dsina — Lcosa (B-4) 


The gravitational forces can be resolved into body axis components such that 


X, = —mgsin@ (B-5) 
Y, = mgsin ¢cos6 (B-6) 
Zg = mg cos ¢cos8 (B-7) 


These equations will be used in the derivations of the V, @, and B equations. Thus, the total forces in the 
body axes can be defined and expanded as 


LX = X;— Dcosat+ Lsina — gmsiné (B-8) 
LY = Yr +Y + gmsin ¢cos6 (B-9) 
“Z = Zr — Dsina — Lcosa + gmcos¢cosé (B-10) 


B.2 Derivation of V Equation 


Beginning with the definition of V in terms of u, v, and w in equation (1-12), 


V =(wW 40? + ww)? 
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the equation for V becomes 


5 4), 4 a a avy r 
Ba = at +0" + w*) | (B-11) 


which after expanding the derivative and cancelling terms, becomes 
. 1 
V= 7p (uit + vb + ww) (B-12) 


By substituting the definitions for u, v, and w from equations (1-9) to (1-11) and cancelling terms, equa- 
tion (B-12) yields : 


V = tcosacosf + tsinG + wsinacos ZG (B-13) 


The definitions for u, #, and w in equation (B-1) are now used with equation (B-13) to give 


Y= Seer ae Xp + Xg) + cosacos A(rv — qw) 


+ me (Yat+ Yr + Y,) + sin B(pw — ru) 


sin a cos G 
+ pa AS TEN ELEY 


(Za + 27 + Ze) + sin acos B(gu — pv) (B-14) 


Expanding (B-14) in terms of equations (B-2) through (B-7) and cancelling yields 


ve “[- Dcos 8 +Y sin 8 + Xp cos acos@ + Yr sini + Zrsinacos A 
— mg(cos a cos Z sin 6 — sin sin ¢ cos 8 — sin a cos 6 cos ¢ cos 4)] 
+ rucosacos @ — qwcosacos § + pwsin 8 — rusin B 
+ qusinacos Z — pusin acos 2 (B-15) 


Equation (B-15) can be simplified by recognizing that the terms involving the vehicle rotational rates are 
identically zero, which becomes obvious after substituting for u, v, and w in these terms. Thus, the final 
equation becomes 


. 1 
V= 7,1 ~ Deos8 + ¥ sin B+ Xp cosacosf + Ypsin 8 + Zr sin acosB 


— mg(cos a cos Z sin @ — sin B sin ¢ cos 6 — sin a cos 8 cos ¢ cos 6)| (B-16) 


B.3 Derivation of & Equation 
The equation for @ can be derived from the definition of a in equation (1-13), 
1w 


a=tan7™ — 
u 


Taking the derivative of a with respect to time, 


y= a = a, -1 0 
Ot ap (B-17) 
then expanding and cancelling terms, the equation becomes 
1 at oo 
= Fp yr ue — tw) (B-18) 
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Substituting the definitions of u and w from equations (1-9) and (1-11) into equation (B-18) gives 


wcosa —usina 
(ee -19 
. V cos (B29) 


Using equation (B-1) to substitute for w and w and equations (B-8) to (B-10) to define the forces, 
equation (B-19) becomes, after rearranging terms, 


a= sey, + Zp cosa — X7 sina + mg(cosa cos ¢ cos + sin asin 8)] 
Vmcos $ 
+ V rice cosa — pvcosa — rvsina + qwsina) (B-20) 


which after substituting for u, v, and w from equations (1-9) to (1-11) and combining terms gives 


a= aS SPED + Zr cosa — Xp sina + mg(cosa cos ¢cos @ + sin asin 6)] 
Vmcos 


+q-—tan 6 (pcosa+rsina) (B-21) 


B.4 Derivation of B Equation 


The equation for # is derived from the definition of 8 as given in equation (1-14), 


epee 1 e. 
B= sin 7 
Taking the derivative of @ with respect to time yields 
> d d . 410 
b= ae =a Vv (B-22) 


which becomes, after expanding the derivative, substituting for V, and cancelling, 
B= J [-itcos asin 8 + bcos 8 — wsin asin f] (B-23) 
Using equation (B-1) to substitute for &, +, and w and equations (B-8) to (B-10) to define the forces, 


; 1 
p= an cosasin § (~Dcosa+ Lsina + Xz — mgsin @) + cos B (Y + Yp + mgsin dcos 8) 
—sinasin 6 (—Dsina — Lcosa + Zp + mg cos ¢cos §)| 


1 
+ yl- cos asin B (rv — qw) + cos 8 (pw — ru) — sin asin B (qu — pv)] (B-24) 
Substituting into equation (B-24) for u, v, and w and rearranging terms yields the final equation 


: 1 - 
B= yp sin8 +¥Ycos#— Xqpcosasin@ + Yycos @ — Zp sin asin B 
+ mg(cos asin @ sin 0 + cos Z sin ¢ cos @ — sin asin f cos ¢cos 8)| 
+psina—rcosa (B-25) 
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APPENDIX C—GENERALIZED DERIVATIVES 


The equations defining the time derivatives of the state variables (derived in sections 1.2.1 to 1.2.4) and those 
defining the observation variables (presented in sections 1.3.1 to 1.3.8) are used to determine the generalized 
partial derivatives of the quantities with respect to a dummy variable €. The purpose of these generalized 
derivatives is primarily to facilitate the derivation of the terms in the linearized equations presented in 
section 2.4; however, these equations have also proved to be useful for computer programs and were used 
to verify the results obtained using LINEAR (see Duke and others, 1987). 


C.1 Generalized Derivatives of the Time Derivatives of State Variables 


Equations (1-39) to (1-41) define the rotational accelerations of the vehicle. These equations are used to 
determine the generalized derivatives of these ee 


aE = detl has BE +h= BE +h GE Day ee OE +I, pe + Be 
Op 
— [2pTezte = IeyI3) - qh 7 LyzIa -D 213) + r(Ieyhy +D yl2 = TyzT3)] = Z 
+ [peel _ Lyzto = D,I3) + 2q( Tye = IeyI3) = (Dr = Ieyle + IpyI3)| a 
é] 
_ [pUTryls + Dytn =; IyzI3) + g(Dh Es Igylo + Ip2I3) + @r(Lyzh a Izz12)] xt (C-1) 
o(4) ane tl { OL OM ON OLty OMr ONT 
Eder b= BE +1,—> BE +15 DE +h OE +h>=> OE + ls ae 
0 
— [2pUa 24 = Ieyts) -— qUezt2 = Lyzl4 —D zs) + age aryl + D yta baz Iyzts)] i 


+ [p(Te2lo _ Lye Tq —D 21s) + 2q(Lyzte _ IeyIs) = r(D; In on Inyl4 + mies ag 


“pda he Dio, 14 Os sie) Fore A 4 (C-2) 


at) 1 {1, aL 8M. AN. Oly | OMy aNy 
é  detI \° 0€ tis oe a€ # ds Oe Poe 6 Oe dg ae oe OE 
_ [2pULnels = IzyJe) = qezI3 = Lyzls = D,I¢) + roy + DyIs = IyzIe)| e 


é 
+ [p(Ie2Ts — Iy2Is = DI) + 2q(IyzTs oo IryIe) = r(DzI3 — Layls + IzzI6)| oF 


té] 
= [pUryl3 + Dyts az IyzIe) + q(D_I3 ee” Jey Is + Trz Ig) o 2r(IyzI3 = Izz15)] zi (C-3) 


The quantities 11, Ip, I3, Is, Is, Is, Dz, Dy, Dz, and det I are defined in equations (1-32) to (1-38) and 
(1-42) to (1-44). 


Equation (1-50) defines the decoupled rotational accelerations of the vehicle (p’, q’, and r’), which are 
used to determine the generalized derivatives of the decoupled quantities: 


Av) _ 1 [OL | abr _ 04 
ae ak Fe (rhe ng payee ge + Plas + 1y + 2ales ~ rls) 3! 
= (pley — ql, + 2rlg, + qiz) =| (C-4) 
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a7) _ 1 dM ap aq 
€ = ra = 3 + FY; = (rie + qlyz + 2plyz = rIz) 3 + (rIey — plyz) ae 
Or 
— (ple — Gay — 2tIz2 — plz) | (C-5) 
air’) _ 1 [ON | ONr 
x = Toe + Be + (qlo + rly, + 2plzy - aly) se zt (ple — rly, — 2qIpy — phy) 5! 
oO 
= (qlez =! plyz) ae (C-6) 
0€ 


Equations (1-58) to (1-60) define the translational accelerations of the vehicle. These equations are used 
to determine the generalized derivatives of these quantities: 
av) 1 


aE =={- cos) Fe + cosa cos s “22 


OYT 


Ding os DE 


Ee chin wedge + sin 6 —— 


OZy 
d€ ag OF 


+[- Xrsinacos f+ 6 ea 


+ cos # cos ¢ cos a cos ayee 


7 
+ [Dsin 6+ Y cos 6 — Xzsin § cosa + Yr cos — Zypsinasin B 
op 


+ mg(sin 6 cos asin 6 + cos @ sin cos B — cos 6 cos d sin asin 8)] BE 


— mg(—cos@cos dsin 8 + cos 0 sin dsin cos 8) 5° 


— mg(cos@ cos acos § + sin @sin dsin B + sin 6 cos Psin a cos 3) sh (C-7) 


O(a) _ 1 OZr OX ap aq Or 
€ = oa (75 pe tesa St FY; ~ sina GE) — tan cosa 5P + of — ten Pele 


OV 


1 
~{apral L 4 Zreosa— Xrsina + mg(cosd cos ¢ cos a + six Bin a)]} Se aE 


ee j e . 
af {z rr ee sina — X7 cosa — mg(cos@ cos ¢sin a — sin O cos a)| 


+ tan (psin a ~ rcosa) \e 


in { tan 8 


Peery gh + Zpcosa — Xpsina 


ap 
dg 


(sin 8 cos & cos a — cos O sin a) z (C-8) 


+ mg(cos@ cos écos a + sin @ sin a)] — —aqe cosa + rsin ay} 


cos # sin ¢ cos a) 


- (rea ie loa 
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—S = — sing 5 hagas ~ cosasins “27 + cos FF sin asin =e 


0g aE 


ov 
+ mg(sin @ cos asin 8 + cos @sin ¢cos 8 — cos Ocos psin asin f)| DE 
+ {ltr sina sin @ — Zr cosasin Z + mg(— sin @ sin asin G — cos @ cos é cos asin ()] 


+ pcosa+ rsin a} 5 ae 


+ —,[D cos 6 ~ ¥ sin 8 — X-rcos acos 6 — Yrsin 8 ~ Zr sin acos 3 


+ mg(sin 8 cos a cos 3 — cos @ sin dsin 8 — cosé cos dsin a cos 3)] o 


+ Z (cos 6 cos ¢ cos 3 + cos @sin @sin asin B) a 


+ Z(cos @cosasin ZG — sin @sin dcos G + sin @cos dsin a sin es (C-9) 


ae 


Equations (1-66) to (1-68) define the vehicle attitude rates. These equations are used to determine the 
generalized derivatives of these quantities: 


4) 9 ae P + sin db tan@ st + cos dtané xt (qcos ¢tan @ — vee ie 

+(qsin dsec” 6 + rcos ¢ sec? 6) - (c-10) 
= cos¢ Fe ding oo - — (qsin ¢ + rcos ¢) a (c-11) 
oe = sin psec d Fe + cos psec i + (qcos dsec @ — rsin dsec 0) a 


+ (qsin dsec # tan @ + r cos ¢sec 6 tan @) a (C-12) 


Equations (1-72) to (1-74) define the earth-relative velocities of the vehicle. These equations are used 
to determine the generalized derivatives of these quantities: 


a) = [cos 6 cos asin 6 — sin B sin d cos 6 ~ cos § sin a: cos ¢ cos 6] e 


— V(cos Z sin asin @ + cos 6 cos a cos ¢ cos #) a 


op 


— V(sin 6 cos asin 6 + cos B sin d cos 6 — sin 4 sin acos ¢cos 6) aE 


~ V(sin 6 cos ¢ cos 8 — cos #sin asin ¢ cos 9) a 
‘ ‘ , . : 00 
+ V(cos B cos acos 6 + sin Bsin dsin @ + cos B sin acos sin 0) DE (C13) 
ah [cos 8 cos a cos 6 cos p + sin B (sin ¢ sin 0 cos — cos @sin p) 


+ cos #sin a (cos dsin @ cos H + sin dsin p)] e 


0a 


0g 


— V{[cos 4 sin a cos 8 cos 4 — cos 8 cos & (cos $ sin 9 cos % + sin dsin ~)] 
~— V[sin 6 cos a cos 6 cos p — cos # sin ¢sin 8 cos p — cos dsin w 
+ sin J sin a (cos fsin 6 cos + sin dsin w)] a 


+ V[sin 6 (cos ¢sin @ cos # + sin d sin ~) — cos f sin a (sin p sin 8 cos  — cos ¢sin ~)] oe 


dg 

é ‘ ; ; 06 

— V[cos 8 cos a sin @ cos  — sin 8 sin ¢ cos 6 cos y — cos B sin a cos ¢ cos 6 cos w] DE 
— V[cos 6 cos a cos O sin y + sin J (sin dsin Asin y + cos @cos yp) 

+ cos 3 sin a (cos ¢sin O sin p — sin ¢ cos p)] ia (C-14) 


ay) 


OE = [cos 8 cos e cos @ sin # + sin B (cos dcos p+ sin ¢ sin 6 sin w) 


+ cos Z sin & (cos dsin @ sin 4 — sin ¢ cos %)] a 


— V[cos 6 sin a cos @ sin w — cos f cosa (cos ¢sin Osin ~ — sin dcos )| i 
— V[sin 6 cos a cos @ sin  — cos 8 (cosécos yp + sin dsin 4 sin ) 
+ sin Z sin a (cos dsin @ sin % — sin d cos p)] ne 
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— V[sin 6 (sin ¢cos p — cos dsin @ sin p) + cos A sin a (sin dsin 4 sin + cos d cos )} i 


06 
— V(cos Z cosa sin #sin % — sin B sin ¢ cos 6 sin y — cos B sin a cos ¢ cos 6 sin p) DE 
+ V[cos @ cos a cos 6 cos p — sin B (cos dsin p — sin sin O cos }) 


+ cos @sin @ (cos dain 0 cos + sin dsin p)] e (C-15) 


C.2 Generalized Derivatives of the Observation Variables 


The vector equation (1-90) defining the body axis kinematic accelerations is used to determine the gener- 
alized derivatives of the individual body axis accelerations: 


Oar) _ col ale aD aL 
Be gem OE COS & —>— DE + sina ae 7 + (Dsina + L cosa) Fe 5 dekciep =" ze (C-16) 
Olay k) _ 1 Se oY 0g j ee 
OE ~ goml| OF tie DE =r + 9mcosé cos ¢ => BE GT NGS Oe (C-17) 
O(azx) tae 1 Es < OD OL . ba 
bf ~ goml d€ sing De cos & De RB NE) a 
—gm cos @ sin } ee gmsin 0 cos ¢ | (C-18) 


Vector equation (1-91) defines the output of body axis accelerometers at the vehicle center of gravity 
and is used to determine the generalized derivatives of the individual body axis accelerometers: 


Alar) _ 1 OXT oD 
Spe gett | Oe. ~ cosa => bE + sino SF z+ (Dsina + Loose) Be -,| (C-19) 
A(ay) _ 1 (= =) : 
BE ~ gom OE * BE oy 
ae = =a TOE sina FP - cosa Fe — (Decosa — Lsina) Fl (C-21) 


Using equation (1-93), the generalized derivative of the output of a normal accelerometer at the vehicle 
center of gravity can be expressed as 


Oe at tel OAT aga OU 
= DE 


BE Gai DE + cosa 5 + (Dcosa — Lsin a) “ (C-22) 


The vector equation (1-95) defining the output of orthogonal accelerometers aligned with the body axes 
but displaced from the vehicle center of gravity is used to determine the generalized derivatives of these 
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quantities: 


Mae) Bae 4 3 28 + ate — Dare) + (pe Orne) 4 ny, 2) (6 
ae = OE i (q¥a + rzz) De t (PU. 2qzz) De t (Pz 2rzz) pe t 22 De Y Be (C-23) 


1 ai NOE ine ae OF nc OR 2 | : 
: (2pyy Wy) Fe (pty + r2y) 5 (qzy 2rtv) Be + 7 BE ty 5g (C-24) 


ag 0 gg 
a2) _ Oaz 1 7 dp 7 oq or ap ai 
Ok 0€ go [cere rz z) O€ + (2gz, TY) 0€ (pz, + Tz) ae Yz ae + Vz =| (C 25) 


Equation (1-96) defines the output of a normal accelerometer aligned with the z body axis but not located 
at vehicle center of gravity, Gp,i. This equation is used to determine the generalized derivative of ans: 


ap dq or 8, Bi 
+ go ez = TZ) ae + (2922 ri TY) dg (pez obi qYyz) FY Ye BE + &;z 5 (C-26) 


O(ani) _ Gan 

0 OE 

In equations (C-20) to (C-23), the partial derivatives of the vehicle rotational rates with respect to the 

dummy variable € are defined by equations (C-1) to (C-3). The partial derivatives of the outputs of the 

body axis accelerometers at the vehicle center of gravity are defined by equations (C-16) to (C-19). In these 

equations, as before, the subscripts x, y, and z refer to the 2, y, and z body axes, respectively, and the 

symbols 2, y, and z refer to x, y, and z body axis locations of the sensors relative to the vehicle center 
of gravity. 


Using cquation (1-97), the generalized derivative of the load factor can be defined as 
O(n) _ 1 OL 


€ — mg OF al 


Equations (1-98) to (1-105) define the air data parameters of interest for this report. These equations 
are used to determine the generalized derivatives of the air data parameters: 


a(a) 0.7p0 aT 
eS O_o C-28 
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ie tae AO C-30 
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Le ge a (C-32) 
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M 


(C-33) 


(C-34) 


(C-35) 


In the preceding equations, the generalized derivative of Mach number appears several times. This term 
can be expanded using equation (C-29). 


The definitions of the flightpath-related parameters are presented in equations (1-106) to (1-108). These 
definitions are used to derive the generalized partial derivatives of the flightpath-related parameters: 


Ay) _ 


ag 


A(fpa) _ 


0g 


ah) 
OE 


1 [hav ak 
~ (V2=A2)U2 | V OE © OE 

1 av 

g OF 

= [-a,,cos¢cos# + a, sin pcos 6] a 


+ [az cos @ + ay, sin dsin # + a,, cos psin 8] a 


+ sind Baz — sin dcosé Day — cos @cos @ was 


0€ 0€ 


(C-36) 


(C-37) 


(C-38) 


The partial derivatives of altitude rate A and velocity rate V that appear on the right-hand side of these 
equations are defined in equations (C-13) and (C-7), respectively. The partial derivatives of the body axis 


accelerations appearing in equation (C-38) are defined in equations (C-16) to (C-18). 
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Using equations (1-109) and (1-110), the generalized derivatives of the energy-related parameters are 
defined. The partial derivatives of altitude rate and velocity rate appearing in equation (C-40) are defined 
in equations (C-13) and (C-10), respectively: 

O(Es) _ V OV , Oh 


dE g OF + 9 (C-39) 


aR) _Vov var, ak 


“ae ~ 9 d€ 1G 8G DE 


(C-40) 


The derivatives of the force parameters, lift (eq. (1-111)) and drag (eq. (1-112)), are defined in sec- 
tion D.1. The generalized derivatives of the normal force (eq. (1-113)) and the axial force (eq. (1-114)) are 
presented in terms of the generalized derivatives of the lift and drag forces: 


O(N) _ OL... aD : Oa 
pe = cose DE + sina aE —(Lsina sel! 73 (C-41) 
oe sina Fe beosa SP ~ (Leosa + Dsina) 2 (C-42) 


The body axis rates are defined in equations (1-9) to (1-11). The time derivatives of these terms are 
defined in equations (1-115) to (1-117). These equations are used to derive the generalized derivatives of 
the body axis rates and accelerations: 


we = cone cos6 Fe ~ Vsin cos 8 22 — V cosasin g 28 (es3) 
= = sin 8 Fe + V cos 0 3 (C-44) 
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(wb) _ 7 = (= - re ee “ae) 
OE DE — sina BE cos a bE V sin 8 22 4. V cosacos $ 5 
+ (qcos a cos § — psin) Se els aie tae BE 


— (qV cosasin 6 + pV cos gy 2 —gcos@sing ae gsin @ cos ge ie (C-48) 
The outputs of various instruments displaced from the vehicle center of gravity are defined in equa- 

tions (1-118) to (1-121). These equations define angle of attack, angle of sideslip, altitude, and altitude rate 

instrument outputs. The generalized derivatives of the quantities are based on these equations: 


O(a:) _ _Ya OP, ta OY _ (Ita = PY) OV | Da 
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: . 8 
+ [- 6(z; sin 6 — y; sin dcos@ — z, cos pcos 8) + H(y;, cos dsin O — cag x 


— (yj, cos cos @ — 2; sin pcos) Sf c+ (xj, cosO + yj, sin dsin 6 + z, cos Psin ee 5 


dh 

=e C-52 
+36 (C-52) 
The generalized derivatives of bank angle rate, pitch attitude rate, and altitude rate with respect to the 
dummy variable £ are defined in equations (C-10), (G-11), and (C-13), respectively. 


The final set of observation variables is defined in equations (1-122) to(1-125). These equations, defining 
total angular momentum and the stability axis rotational rates, are used to determine the generalized 
derivatives of these quantities: 


MOD (Lap — Foyt ~ Fost) FE + (Iya — Tey — Ir) $2 + (Let = Tos In) $e (C53) 
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APPENDIX D—EVALUATION OF DERIVATIVES 


The generalized partial derivatives presented in equations (C-1) to (C-56) contain partial derivatives of the 
state variables, thrust forces, and total aerodynamic forces and moments with respect to the dummy variable 
é. In this appendix, these partial derivatives are defined with respect to specific state, time derivatives of 
state, and control variables. The derivatives of atmospheric parameters are also discussed. 


D.1_ Preliminary Evaluation 


First, the partial derivatives of the state variables with respect to the state, time derivatives of state, and 
control variables are considered. All partial derivatives of the state variables with respect to the state 
variables are either equal to zero or unity. Thus, 


dp_0¢_or_ OV _d4_98_ 0 8 Ob _dh_% Ww, my 
Op O¢ Or AV da Of O6 G6 Ab Gh Ax dy 
and all other derivatives of state variables with respect to state variables are equal to zero. The partial 
derivatives of the state variables with respect to the time derivatives of the state variables (@ and B, in 
particular) are equal to zero. This is also true of the partial derivatives of the state variables with respect 


to the control variables. 


Second, the partial derivatives of the aerodynamic forces and moments with respect to the state, time 
derivatives of state, and control variables are evaluated. Using the definitions of the force and moment 
coefficients presented in appendix A, the partial derivatives can be explicitly evaluated in terms of the 
stability and control derivatives. 


D.1.1 Rolling moment derivatives.— 


= . ae (D-2) 
a4 vt 09 
= Ce (D-4) 
ao = SbpVCr + qS5bC2, (D-5) 
oh = §5bCr, (D-6) 
= = G5bCe, (D-7) 
a = 550V°Cr as LiGSie,: (D-8) 
3 . ae (D-9) 
oi = ee, (D-10) 
ve = gSbC%,. (D-11) 
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D.1.2 Pitching moment derivatives.— 


= StpVCm + GStCmy (D-15) 
os = Cra (D-16) 
Oy = GST, 7) 
ae = 55070 oe + GStCm, (D-18) 
ae 2 IE Cm (D-19) 
4 = Cn, (D-20) 
= 1S, (D-21) 
D.1.3 Yawing moment derivatives.— 
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oe Ou (D-29) 
ay = Gy (D-30) 
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D.1.4 Drag force derivatives — 


D.1.5 Sideforce derivatives.— 
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D.1.6 Lift force derivatives.— 
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Next, the partial derivatives of the powerplant-induced forces and moments with respect to the state, 
time derivative of state, and control variables are considered. The partial derivatives of the powerplant- 
induced forces and moments are assumed to be zero except for moments taken with respect to the body 
axis rates (p, q, 7), moments and forces taken with respect to the velocity and velocity orientation terms 
(V, a, @), and forces taken with respect to the control variables. These terms, assumed to be nonzero, are 
taken as primitives and not evaluated further. Thus, using F, to represent a powerplant-induced force (X7, 
Yr, and Zp) and M, to represent a powerplant-induced moment (7, Mr, and Nr), 


OFp _ OFp _ OFp _ OF _ OF _ OFp _ OF _ OF _ OF _ (D-62) 
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ad 00 Ow Oh Ox Oy 08; 


d 
a OF, ON) OF, OF, 2M, OM, OM, OM, OM, OM, 


OV’ Oa’ OB’ 06; Op’ dg’ dr’ OV’ da’ OB 
are taken as primitives and not evaluated further. 


The final set of partial derivatives to be discussed are the derivatives of atmospheric parameters with 
respect to the state, time derivative of state, and control variables. In this report, all atmospheric parameters 
are assumed to be functions of altitude only. Thus, except for 


OT Op Op Ova 

Oh? dk? Oh? 4 Bh 
all derivatives of ambient temperature, density, viscosity, and ambient pressure are assumed to be equal 
to zero. The nonzero quantities listed previously are dependent on an atmospheric model. Clancy (1975), 
Dommasch and others (1967), Etkin (1972), and Gracey (1980) present discussions of atmospheric models. 


In this report, the quantities will be taken as primitives and not evaluated further. 
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D.2. Evaluation of the Derivatives of the Time Derivatives 
of the State Variables 


The generalized derivatives of the time derivatives of the state variables are defined in appendix C, equa- 
tions (C-1) to (C-15). In this section, these generalized derivatives are evaluated in terms of the stability 
and control derivatives, primative terms, and the state, time derivative of state, and control variables. In 
this section, the notation 0(x;)/0z; is used to represent the more correct notation 0f;/Ox; that is employed 
in the discussion at the beginning of section 3. This notation is used because there is no convenient no- 
tation available to express these quantities clearly—particularly not the usual notation employed in flight 
mechanics texts such as Etkin (1972) and McRuer and others (1973). The notation that defines quantities 
such as Ly = O(p)/Op and M, = O(¢)/0q is misleading in this context because the definitions of those terms 
(such as L,, M,) are based on assumptions of symmetric mass distributions, symmetric aerodynamics, and 
straight and level flight, and additionally do not include derivatives with respect to atmospheric quantities. 
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D.2.2 Pitch acceleration derivatives.— 
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D.2.3 Yaw acceleration derivatives. 
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D.2.4 Decoupled roll acceleration derivatives.— 


£ievbot B06 = 1)| 


= Ley Do + go(ly re I,) Fe 2Iyer0 


O(p’) _ 1 /gse? OL 

dp Ts Ea Oe “Fp ~ favt0 + Teeao| 
ap’) 1 EE aly 

“8q. oT L BV, Cat Og 

aw) _ 1 [2s OLy 

Or ‘I, | 2Vo Co + Ge 

av) _ 1 . OL 
Or L Sb(pVoCe + GCu,) + OC 
ar) 1 (« i) 

aa. Te Be 

ap!) _ ale if 2 OLr 

“Be 7 1, (e+ ap) 

ap’) 

oan 

Ap!) 

"9g 

ap!) _ 

ap 

Ap’) Sb op 

Be = 7, (g¥BCe 5p + aCe) 
d(p') 

oc u 

Ap!) 

oy : 

Ap!) _ GSbe 

Oa -Wol, 

a(p') _ gS? 

0B Vol, 

A(p') _ gb 

66; ~ T, Coy 


(D-106) 
(D-107) 


(D-108) 


(D-109) 
(D-110) 
(D-111) 
(D-112) 
(D-113) 
(D-114) 
(D-115) 
(D-116) 
(D-117) 
(D-118) 
(D-119) 
(D-120) 
(D-121) 
(D-122) 


(D-123) 


D.2.5 Decoupled pitch acceleration derivatives.— 


ao = z Ea mp + — 2Iz2Po — Lyz9o + To(Le — I.) (D-124) 
a 1 (Been, 4 PM taro It) (D-125) 
AO = r Cn, + oat + po(Iz — Ie) + Toyo + 21.70 (D-126) 
aa) = a [$a(e¥eCm + GCmy) + | (D-127) 
20 He om 
“ee =0 (D-130) 
a) = (D-131) 
= a, (D-132) 
aa) =0 (D-134) 
oa =0 (D-135) 
ow] 2 
a I IRE a (D-136) 
ne 
4 7 arom (D-137) 
af) = One (D-138) 
D.2.6 Decoupled yaw acceleration derivatives. — 
“ = : Ee + ee + 2Inypo + qo(le — Ty) + her (D-138) 
oo = z [Seen + oe + poz — Ty) ~ 2Lzy%o — Testo] (D-140) 
“4 (F ae a shen 10) (D-141) 
AE) = F [s6(0V0Cn + ACmy) + (D-142) 
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Ja 1, (15m + FE) (D-143) 
oT . z (as bCng + oP) (D-144) 
ae aM (D-145) 
oP au (D-146) 
i me (D-147) 
“Wh Te (Zn 5p + 40m) (0-48) 
=! (D-149) 
o =e (D-150) 
7 - eC (D-151) 
oF = ren, (D-152) 
x _ FC (D-153) 


D.2.7 Total vehicle acceleration derivatives — 


AV) _ 456 


Bp = BVam (7 cos Bo Cp, + sin Bo Cy,) (D-154) 
av qSt . 
ore = eae cos Bp Cp, + sin Bo Cy,) (D-155) 
av gSb . 
ay) = cae cos Jo Cr, + sin Bp Cy,.) (D-156) 
av) 1 _ : 2 
Ay = al S'cos Bo oe + q@Cp,) + Ssin ‘ Sieg) + GCy,) 

oY; 
+ cos ag cos pee a ee sin Qg cos pees a ise sin Bp | (D-157) 
av 1 
an) =— =| — FS cos Bo Cp, + G5 sin Bp Cy., + cos ao cos wo 


O27 OY> 
+ sin ao cos Bp —— Fa sn eRe auens 


+ Zp cos ag cos fy + mg(sin & sin &g cos By + COs Oy cos dg COs ap COS 6o)| (D-158)} 


1 
— = |a5(- cos Jo Cp, + sin 89 Cp + sin Bo ae + cos Bo 


a Blt 4. Go stg 
COS Ao COS Vag — dB sil. Mo cos (areprae 55 sin O75. oe 


— Xp sin Bp cos ag ~ Zp sin ag sin Bo + Yt cos fo 


+ mg(sin 99 cos ag sin Bo + cos Op sin do cos Bo — cos Oo cos go sin ap sin Gq) 


g(cos @ cos ¢o sin Bp — cos Ag sin $o sin a cos fo) 


g(— cos Ao cos ao cos Zo — sin Mo sin ¢o sin Go — sin Mo cos do sin an cos fio) 


0 
=| bo (5V8C> 5h IC, ) Bo (Svecy C 
he eee 0 Dap t4 D, } + sin Bo 0 Vag e vn) 
0 
0 
SE . 
= Wan cos Jo Cp, + sin Bo Cy, ) 


GS 3 
som =o cos By Cp, + sin Bo Cy;) 


Ss 
“(- cos Bo CDs, + sin Bo Cy,,) 


1 OZy 
+— = (cos 2905 boa + sin Bo oe + sin a cos Bo D5; oe 


D.2.8 Angle-of-attack rate derivatives — 


~ 2V2m cos Bo 


~ 2V2m cos Bo 


g5b 
qs CL, — tan Bo cos ag 


ae eee 


____45b . 
2V2m cos fo Pip = tan foein Ge 

ae OZy . OXT 

mVi cone {s(WopC + GCL) — cosa WV + sin ag ae 


1 
+ GIS Cu + Zpcosag — XT sin ag 
0 


+ mg(cos 99 cos ¢g Cos Ag + sin Op sin ag) 


(D-159) 
(D-160) 
(D-161) 
(D-162) 
(D-163) 
(D-164) 
(D-165) 
(D-166) 


(D-167) 


(D-168) 


(D-169) 
(D-170) 


(D-171) 


(D-172) 
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8(a) 1 OXy 


lsc, — COS Qo oer + sin ag + Zy sin ag + Xp cos ao 


a ne mVo cos Jy da 0a 
+ mg(cos 9 cos ¢o sin ao — sin My cos ag) 
+ tan Bp (po sin ag — ro cos ao) (D-173) 
O(a) _ 1 {a 7 OZ, i, Oxy 
ym = es GSC, — cos ag OB + sin ao yn 
— tan Bo [-GSCL + Zz cos ag — Xz sin ao 
+ mg(cos 69 cos }p cos ag + sin Op sin ao)]} 
- aa Po COS 9 + Tp Sin Ag) (D-174) 
0 
O(a) g : 
ON eh D-175 
a6 Ween cos 49 sin $p cos ao ( ) 
O(a) _ g . _ F : 
ae “Vocorag en 89 cos $p cos a — cos Gp sin ao) (D-176) 
~ 90 (D-177) 
O(a) _ Ss lian Op, . ) : 
ao) =0 (D-179) 
He) _ (o-a 
y 
O(a) _ gGSé 
da ~~ -V2mcos Bo Cha Date) 
O(a) _ gSb 
ap a ~2V2m cos Bo La (D-182) 
a) : aZp or) . 
06; = mVo cos By [-75C1, + COs A On. — SIN Qo Os, (D-183) 
D.2.9 Angle-of-sideslip rate derivatives — 
a6 GSb .. : 
2 = sre (sin Bo Cp, + cos By Cy,) + sin a (D-184) 
0 
a(6 gse 
oe = spa (sin Bo Cr, + cos By Cy,) (D-185) 
0 
a(8) = 285) ie Bo Cp, + cos Bo Cy,) — cos ag (D-186) 


Or 2Vem 


a(6 1 
a avs [s sin fo (pVo + ¢Cp,,) + S cos Bo naa + gCy,) 


OV 
Oxy OZ 
— ¢€08 ao sin Sp —— a + cos Jo oo — sin ag sin fo a 
1 
Ve [GS (sin Bo Cr + cos Bop Cy) — Xz cos ag sin Bo 
0 
+ i cos 89 — Zr sin ag sin Bo] (D-187) 
ate 
of) = AV —ae |Bs(ein Bo Cp, + cos Bo Cy,,) — cos ag sin fo % T 4 cos $ by ST 
— sin a sin Bp ofr + Xp sin Qp sin Bo — Zr COs Ao sin Bo 
— mg(sin 4 sin ag sin Bg + cos A cos do cos ag sin 6o)| 
+ po cosag + rosin ao (D-188) 
OEY Ne, qS [si C C Cot+C 
BG = VES sin Bo (Cp, — Cy) + ees p+ Cy,)} 
— cos ao sin Bo oaE +c 0s Bo oe — sin a sin fo ee 
— XJ cos ag cos Bo — Yz sin Bo — Zr sin ag cos Go 
+ mg(sin 69 cos ag cos Bp — cos 9 sin do sin Bo — cos Op cos Pp sin ao cos bo) (D-189) 
a = 2 (cos 9 cos do cos Bg + cos 89 sin Po Sin ao sin Bo) (D-190) 
0 
“f) = A cos 9 cos ag sin Bo — sin 9g sin go cos Bo 
0 
+ sin p cos do sin ag sin fo) (D-191) 
ats 
ee =0 (D-192) 
1) en 1 dp 1 dp | 
Ph = mVo [sin Bo (Greco 3 oh + 70D, ) + cos Bo (Svecy dh + cv, (D-193) 
ate 
8) = (D-194) 
AtL 
“) 0 (D-195) 
Bie 
“) = ign aoe (i 8o Cp, + cos Bo Cy,,) (D-196) 
ve 
aD = = ie —>—(sin Bo Co, + cos Bg Cy; ) (D-197) 
AB) _ 1 
ae) = Vp [25 (sin Fo CDs, + cos fio CY, ) — cos ao sin fo St et cos fo 5 — 
4 
: P OZy 
— sin ag sin fo a} | (D-198) 
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D.2.10 Roll attitude rate derivatives — 


—— = sin do tan 9 


(4) 
Or 


a(¢) 


d¢ 


Z o|— 
Q- Sle 
SS S 


= C08 dg tan Ap 


= qo cos dp tan Op — ro sin dp tan A 


= go sin $9 sec” % + To cos do sec? Io 


D.2.11 Pitch attitude rate derivatives— 


08) _ 
op >" 
— = cos Go 


0) = —sin do 


OV 


(D-199) 
(D-200) 
(D-201) 
(D-202) 
(D-203) 
(D-204) 
(D-205) 
(D-206) 
(D-207) 
(D-208) 
(D-209) 
(D-210) 
(D-211) 
(D-212) 


(D-213) 


(D-214) 
(D-215) 
(D-216) 


(D-217) 


a4) _ 4 (D-218) 


da 
2) ai) (D-219) 
1 = —go sin do — Tro c0s $o (D-220) 
a) 49 (D-221) 
a x, (D-222) 
uU) =0 (D-223) 
a -0 (D-224) 
4) = (D-225) 
a6) = (D-226) 
— =0 (D-227) 
se) =, (D-228) 
D.2.12 Heading rate derivatives.— 
a) ae (D-229) 
A) = sin go see (D-230) 
2) — cos 500 Ey (D-231) 
5) 2% (D-232) 
oe) a (D-233) 
uy a5 (D-234) 
a) = qo cos Po sec Oo ~ ro sin do sec A (D-235) 
ae) = qo sin do sec O tan 9 + To cos do sec Oo tan A (D-236) 
a) sn (D-237) 
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A) = 0 (D-238) 
0 (D-239) 
ta =" (D-240) 
8) = 9 (D-241) 
Oi. (D-242) 


a 


D.2.13 Altitude rate derivatives.— 


ae a (D-243) 
2 ~ (D-244) 
eo ve (D-245) 
zs = sin 9 cos fig cos Qo — sin do cos Ay sin By — cos bg cos Op cos fo sin ag (D-246) 
th) = —Vo(cos Ap sin ag sin 8 + cos Bo cos ag cos do cos Bo) (D-247) 
) = —Vo(sin Ao cos ag sin 89 + cos fo sin }o cos 9 — sin fo sin ag cos $y cos 99) (D-248) 
au) = —Vo(sin Bo cos $9 cos Oy — cos Bo sin ag sin do cos 89) (D-249) 
at = Vo(cos Bo cos ag cos 8 + sin fo sin go sin 8 + cos Ao sin exo cos do sin om) (D-250) 
se 7" (D-251) 
- ai (D-252) 
oe a (D-253) 
o = (D-254) 
We = (D-255) 
% =e (D-256) 
m= (D-257) 


D.2.14 North acceleration derivatives.— 


D.2.15 


=0 
=0 
=0 


= cos fg COs Qo Cos Op cos Ho + sin Bo (sin do sin % cos Yo -- cos go sin Ho) 


+ cos Bg sin ag (cos do sin Og cos Hq + sin Po sin 9) 


= Vo[ cos Bo cos ao (cos $9 sin A cos Ho + sin do sin Yo) — cos Bo sin ao cos Oo cos to] 


= Vo| cos fo (sin do sin 89 cos yo — cos hg sin yo) — sin Bp COs ag COS Ap COs Yo 


— sin Bo sin ag (cos $o sin Op cos wo + sin do sin Ho)] 
Vo[ sin So (cos go sin 9 cos Ho + sin do sin Yo) 

+ cos J sin a (cos do sin Ho — sin Po sin Op cos %o)} 
Vo( sin Go sin ¢o cos 89 cos Ho — cos Bo cos ao sin 8 cos Ho 


+ cos fg sin ag cos dg cos Jo cos ¥p) 


= Vo|[ — cos Bp cos ag cos 9p sin Ho — sin fp (sin go sin Op sin Wp + cos og cos Yo) 


— cos Bo sin a (cos o sin Op sin #o — sin dq cos Ho)] 


(D-258) 
(D-259) 


(D-260) 


(D-261) 
(D-262) 


(D-263) 


(D-264) 


(D-265) 


(D-266) 
(D-267) 


(D-268) 
(D-269) 
(D-270) 
(D-271) 


(D-272) 


(D-273) 
(D-274) 


(D-275) 
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= = cos sin to cos G9 cos ag + sin Bo (cos do cos Yo + sin do sin Mo sin yu) 


OV 
+ cos fig sin ag (cos do sin 89 sin g — sin dg cos thy) (D-276) 
oa) = Vo[cos Bo cos a (cos do sin 8p sin Yo — sin Pp cos Yo) 
— cos Jo sin a COs Ig sin Yo} (D-277) 
ae = V| cos fo (cos go cos Ho + sin do sin 8p sin Wo) — sin Jo cos ag cos 9 sin Wo 
— sin Bo sin ao (cos go sin % sin Ho — sin gg cos Yo)] (D-278) 
3 = Vo[sin fo (cos do sin 8 sin Wo — sin do cos Yo) 
— cos Bo sin a (sin go sin Hp sin Wo + cos $o cos yo)] (D-279) 
aD) = V( sin Bo sin dp cos Op sin Wp — cos Bq cos ap Sin Op sin wo 
+ cos fo sin ag cos Pp cos Op sin to) (D-280) 
a = Vol cos Bo cos ag cos Op cos Yo — sin Bo (cos go sin Yo — sin ¢o sin Op cos Vo) 
+ cos Jo sin ap (cos $0 sin 4 cos Yo + sin do sin Yo)] (D-281) 
ay) 
= D-282 
Oh : ( ) 
ay) 
pala SE Are -283 
ae = 0 (D-283) 
ay) 
pS ee D-284 
Be = 0 (D-284) 
Oy) _ 
a =O (D-285) 
9) _ 4 (D-286) 
dg 
ay) _ 
06, = 0 (D-287) 


D.3 Evaluation of the Derivatives of the Observation Variables 
The generalized derivatives of the observation variables are defined in appendix C, in equations (C-16) to 


(C-56). In this section, these generalized derivatives are evaluated in terms of the stability and control 
derivatives, primative terms, and the state, time derivative of state, and control variables. 


D.3.1 Longitudinal kinematic acceleration derivatives.— 


O(a, qSb ‘ 
~ Wagar cosa Ch, + sin ao Cr,,) oe) 
ee = TG cos ag Cp, + sin ao CL, ) ay) 
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ae = aa - S cos aq (pVoCp + GCp,) + S sin ao (pVoCL + GCL) + a (D-291) 
“Ge = a {as [— cos ao (Cp, — CL) + sin ao (CL. + Cp)] + “ar (D-292) 
a = aa [as(- cos ao Cp, + sin ag CL,) + oa (D-293) 
aa a (D-294) 
ee = —F,c084 (D-295) 
oe =e (D-296) 
Mer) - a — COS Ag (Fsvecns6 + 18Cp,) + sin ao (ssvec.se + SC, )| (D-297) 
ee = (D-298) 
ee) 7 (D-299) 
Sex) = Fen (- cos ao Cp, + sin ao CL) (D-300) 
“ee Ee Feat cos a9 Cp, + sin oo CLs) (D-301) 
aon) > aan —— [as(- cos ao Crp,, + sin a CL,,) + Pal (D-302) 
D.3.2 Lateral kinematic acceleration derivatives.— 
one = Fieger (D-303) 
a Ta (D-304) 
al = Man (D-305) 
“Cu) 7 aan —— (Sp¥eCy + 48Cy, + +r ) (D-306) 
ae jam (180%. + Fe) (D-307) 
ou) == a (a SCy, + ony fea 
wee = a cos 69 cos go (D-309) 
Meu) = = sin 09 sin do (D-310) 
ape =e (D-311) 
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Aout) _ = (5 = SVECy oe iSCy,) (D-312) 
ee %4 (D-313) 
Cee = (D-314) 
LE = 5 icv, (D-315) 
Z = 5 icv, (D-316) 


D.3.3 Z-body axis kinematic acceleration derivatives — 


O(az%) _ gSb ‘. 
a —BVogamn oi ao Cp, + cos a CL,) (D-318) 


(az) qst_,. 
“a = — Bogan ao Cp, + cosa Cy) (D-319) 


Odsk) _ __ ab 
Or ——- 2Vogom 


O(a, 1 
Sok) = “tom S sin a9 (pVoCp + GCp,,) + S cos ag (pVoCL + FCL,) — 
1 dZr 


O(az%) = . 
ae ae = - sin ao (Co, = CL) + cos Qo (CLa + Cp)| + Jom Da (D-322) 
1 O27 


O(a) 
gom Of 


dp 
= sr cos 4 sin ¢o (D-324) 
0 


(sin ap Crp, + cos ag Cy) (D-320) 
a 


(D-321) 


qs 
= -aasin ao Cp, + cos a9 Cig) + — (D-323) 


O(@zk) 
dg 
O(azx) 
06 
O(az) 
OY 


O(az if. 1 Op Op 
XC) = — =| sin Qo (SsveCp 3 Dh + SCD,) + cos a (5 SVECL = ah +a5C1,)| (D-327) 


= ain 09 cos by (D-325) 
Jo 


=0 (D-326) 


Aes) =0 (D-328) 
O(az%) a 
aa (D-329) 


O(a, gSc 

“een = OVae mac ao Cp, + cos ag CL, ) (D-330) 

O(azk) _ _ 55h 
dB  — — 2Vogam 

O(azk) 
06; 


(sin a9 Cp, + cos ao CL) (D-331) 


oz) 


oF (D-332) 


1 ¥ "y 
= ra g5(sin ag CD, + COS Qg Ch,,) - 


D.3.4 x body axis accelerometer output derivatives.— 


oe = Sent cos do Crp, + sin a CL,) 

Pee) ; - : (— cos ao Cp, + sin ao Cr) 

Mae) = Scent cos ag Cp, + sin a CL.) 

iS) — [- S cosag (pVoCp + Coy) + S sin ao (pVoCL + GCLy) + aes a | 
a = = {as[- cos a (Cp, ~ CL) + sin ao (CL, + Cp)] + =— 

Mes) —- ae cos a9 Cr, + sin ao Cig) + aa 

wee! 0 

a 

a 

Kos) os — | — cos ag Gage se fs 15Cp,) + sin a9 (Gsvec. oe + SCL, )| 
a 

es) as ae (—cosag Cp, + sina CL) 

e = sal COS Qo Cp, + sin ag Cy,) 


D.3.5 y body axis accelerometer output derivatives — 


Olay) gb 
Op a 2Vogom Xp 
Alay) _ _7S7_o 
Oq ~ 2Vogom Mg 
O(ay) g5b 
“Or ee ae 
(ay) 


7 avy 
—— (Spvocy a ore +5) 


(D-333) 
(D-334) 
(D-335) 
(D-336) 
(D-337) 
(D-338) 
(D-339) 
(D-340) 
(D-341) 
(D-342) 
(D-343) 
(D-344) 
(D-345) 
(D-346) 


(D-347) 


(D-348) 
(D-349) 
(D-350) 


(D-351) 
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DOB ay cd (acy. + at) 


da gom 0 
A(ay) 1 (a ain) 
ey) te rT 

OB = gom \PC¥e + Gp 
Ody) _ 

age 
O(ay) 

Der ee 
O(ay) -_ 

Op a 
Ody) 1 Op 
Oh es gai am (5 SVS Cy ae F + 45Cy,) 
O(ay) 

4 
Ady) -_ 

Oy M 
O(ay) GSE 

Oa ~ 2Vogom Ya 
O(ay) _ — g&b 

a6 — 2Vogom~ *8 
O(a 1 _ Oy: 

Be ium (290%, + Fe) 


D.3.6 z body axis accelerometer output derivatives.— 


es) 


G5b 

Sa (sin ap Cp, + cosa CL, ) 

qSé 
ers a9 Crp, + cos a CL,) 

GSb 
Wear (sin Oo Cp, + cos a Ch, ) 

1 
aa [s sin ag (pVoCp + GCp,,) + S cos ao (pVoCL + GCL,) - =| 

: az 
{astsin ao (Cp, — CL) + cos a (CL, + Cp)] — =} 

gom 

wis O27 
Fes — |75(sin a0 Cp, + cos a9 CL) — a 
—— {si (sve op + GSC; ) cos € SV2C; Oe reg GSC )] 

pe 1 AQ 3° "0 DSF qolp, | + 0 0 L On qOUL, 


(D-352) 
(D-353) 
(D-354) 
(D-355) 
(D-356) 
(D-357) 
(D-358) 
(D-359) 
(D-360) 
(D-361) 


(D-362) 


(D-363) 
(D-364) 
(D-365) 
(D-366) 
(D-367) 
(D-368) 
(D-369) 
(D-370) 
(D-371) 


(D-372) 


=e | (D-373) 


Ox 
al =a (D-374) 
Ge) = - (sinag Cp, +cosaoCL,) (D-375) 
ee) = -5 a (sin ao Co, + COS &g CL, ) (D-376) 
a = -— [as¢sin ao CD,, + cos ao Cu = a (D-377) 
D.3.7 Normal accelerometer output derivatives— 
aoe = Tcon (sin ap Crp, + cos a CL,) (D-378) 
ae = _. (sin ag Cp, + cosa CL,) (D-379) 
Aen) = Socom (sin ag Cp, + cos a CL, ) (D-380) 
eo) 7 —- [sn a (SpVoCp + FSCpy) + cos ap (SpVoCi + GSCLy) — oer (D-381) 
(Gn) > ai { a8{sin ao (Cp, — Ci) + cos a0 (CL, + Cp)] — ao} (D-382) 
oe = = |as'sin ao Cp, + cos a0 Cig) — a (D-383) 
_ a (D-384) 
ae =" (D-385) 
woe ae (D-386) 
Men) = — [sin a (5sv3cp se a 75Cp,) + cos an (5 sv2c, 2 ve +95C.s)| (D-387) 
ae aa (D-388) 
“Eel =a (D-389) 
en) = (sin ao Cp, + cosa CL) (D-390): 
aa = Teen (sin a9 Cp, + cos ao Cr,) (D-391) 
eel = a G5(sin ap Cp,, + 608 do CLs.) — oe (D-392) 
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D.3.8 Derivatives of x body axis accelerometer output not at the vehicle center 
of gravity.— 


Ades) _ _ qb 


ap BVogom (00820 Cb, + sin ao Cry) + ~ (dot + Tozr) (D-393) 
“Gail = tom cos ao Cp, + sin ag Ci, ) a ~(Pove - 2q022) (D-394) 
*Goii = Seon sO (~ 608 ao Cp, + sin ag CL, ) + ~ (Poze — Wrox) (D-395) 
*eei) = a - S'cos ag (pVoCp + FC) + S sin ao (pVoCL + GCL, ) Paes} — (D-396) 
Meni) = an {ast- cosaq (Cp, — CL) + sin a9 (CL, + Cp)] + S| (D-397) 
"Gesl ern — |a5(- cosag Co, + sin ao Cig) + a (D-398) 
— = (D-399) 
ee a (D-400) 
ee = (D-401) 
oe 7 aa |- COB Go (5 SV8Cp 92 Dh e+ 25Cp, ) +sin ao (Sse. oe + 15C..)] (D-402) 
ru = (D-403) 
ie vi (D-404) 
cs aes (D-405) 
a ~ z (D-406) 
as ~ oa (D-407) 
Aes) - a- cos a Cp, + sin a CL) (D-408) 
aie = Se l- COS Ag Cp, + sin ag Cis) (D-409) 
Me! z an [- G5(~ cosa Cp,, + sin ao Clade Sel (D-410) 


D.3.9 Derivatives of y body axis accelerometer output not at vehicle center 
of gravity.— 


1 
Cy, - Foe oy — Gry) (D-411) 


Maa = ae + ~(Poty + rozy) | (D-412) 
Gui) = Hence + (902 — 2royy) (D-413) 
Novi) s —— (Spvocy +GSCy, + =) (D-414) 
Meus) a (ascv. + a) (D-418) 
oes =9 (D-417) 
Moya) = (D-418) 
ogee -9 (D-419) 
Meus) — (Fsvecy es 75C¥, ) (D-420) 
ous) a4 (D-421) 
wee 0 (D-422) 
Mae) _ -% (D-423) 
oes =0 (D-424) 
Kens) = -* (D-425) 
“Aes = es : (D-427) 
“ul e aaa (ascv,, + Se) (D-428) 


D.3.10 Derivatives of z body axis accelerometer output not at vehicle center 


of gravity— 
(az) gsb_ 1 
a fas —3¥egam ie ao Cp, + cos ag CL) — qe PO — T0fz) (D-429) 
O(az,i) GSE : 1 
Sr an ao Cp, + cos ag CL,) ~ 90 — To¥z) (D-430) 
Olazi) _ gsb 1 
aS “aon ao Cp, + cos ag Ci,) + 9g (Pot? + qoyz) (D-431) 
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D.3.11 


~ gom 
as OZ> 
~~ gom [75(sin 00 Cd, + c08 00 Chg) — Be | 
=0 
=0 
=0 
‘ae dp ee 

= Fa [tm a0 (58VEC FP + a8Cp,) + cosa (IsvEC. 2 + asc, 
= 0 
= 0 

Ya 

go 

_ oz 

go 
=0 
= eae 517-7 (sin ao Cp, + cos ag CL, ) 
gGsb 
= — Fao (tin 20 Co, + COS aig CL, ) 
= : OZ 
see [7s(sin Qo CD,, + cosag CLs.) — 06; 


— |s sin ao (pVoCp + GCp,,) + S cos ao (pVoCy + gCL,) —- 


aa {asilsin ao (Cp, — CL) + cos. ag (CL, + Cp)] — aa 


OZ 
aV 


aa 


7 | 


Derivatives of normal accelerometer output not at vehicle center 
of gravity.— 
a ani GSb 
oe ) = ee (sin ag Cp, + cos ag CL,) er ~(2p022 — Toxz) 
O(Gn,3) qGSé 
ar ve (sin a Cp, + cosa Cy,) + te — Toyz) 
0 ani qb ‘ 1 
(on, ) = Wages ag Cp, + cos ag CL.) — go Pot + qo¥z) 
O(an,i) 
ay mn ——|s sin ao (PVoCp + FCp,,) + S cos ag (pVoCy + GCL) - 
re} ani _ ‘ az 
ion, ) = a { a8{sin ao (Cp, — Ci) + cos ao (Cy, + Cp)] — a 


ar 


(D-432) 
(D-433) 
(D-434) 
(D-435) 
(D-436) 
(D-437) 
(D-438) 
(D-439) 
(D-440) 
(D-441) 
(D-442) 
(D-443) 
(D-444) 
(D-445) 


(D-446) 


(D-447) 
(D-448) 


(D-449) 


- (D-450) 


(D-451) 


; OZ 
nee = aa [7s(sin ao Cp, + cos ap Cig) — a (D-452) 
Aan i) _ D-453 
oe (D-453) 
Aan) _ D-454 
a7 ( ) 
O(an, i) 
Dae D-455 
ce (D-465) 
O(ani) _ 1 . ee) op = ) (5 2 op — )] 4 
Oh Gare [sin ao (55% Cp Dh + qGSCp, | + cos ao 55% CL Dh + gGSCy,, (D-456) 
Alans) _ 0 (D-457) 
Ox 
O(ani) _ 0 (D-458) 
dy 
Alani) _ _ Ye (D-459) 
dp Jo 
Alani) _ te (D-460) 
ag Jo 
Pani) _ 9 (D-461) 
Or 
se) - Ho (sin a9 Cry + .c08 ap C,) (D-462) 
ao = eon (sin ao Cr, + C08 a CL) (D-463) 
nyt = : OZ 
“Gni} = aii [as¢ein ag CrE,, + COs a9 Cis, ) - ae (D-464) 
D.3.12 Load factor derivatives.— 
O(n) _ gS : 
ap BVogme eae) 
O(n) qSe 
canes Aon -4 
Bq DVogme (D-466) 
An) _ — g5b ; 
or 3Vogm br Dae) 
O(n) 1 > 
BV = al SeV0Cr + 28CLy) (D-468) 
Bin), 2/08: : 
ae gmc (D-469) 
O(n) — gS 
>—_ => D-470 
a6 pers ( ) 
An) _ ’ 
30 = (D-471) 
O(n) _ ; 
ao = 9 (D-472) 
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op 9 
On) 1 flag, Op _ ) 
We oe (55¥8 Cu ge + AS8C1, 
a(n) 
oz 
O(n) 
dy 
O(n) Se 
Oa ~ 2Vogm la 
an) gb 
QB  2Wegm ve 
O(n) gS 
05; gmt 
D.3.13 Speed of sound derivatives — 
O(a) _ 
dp q 
Ala) _ 
Oq 
O(a) 
or 
O(a) 
av? 
(a) _ 
da 
O(a) _ 
op 
(a) 
fy ae 
O(a) _ 
oo 
Aa) _ 
ay 79 
O(a) _ 0.770 oT 
Oh ~ poTo(1-4 po/poTo)'/? Oh 
O(a) 
or 
O(a) _ 
ayo 
O(a) _ 
pa 
Aa) _ 9 
B 


(D-473) 
(D-474) 
(D-475) 
(D-476) 
(D-477) 
(D-478) 


(D-479) 


(D-480) 
(D-481) 
(D-482) 
(D-483) 
(D-484) 
(D-485) 
(D-486) 
(D-487) 
(D-488) 
(D-489) 
(D-490) 
(D-491) 
(D-492) 


(D-493) 


Aa) _ 9 (D-494) 


oO 6 (D-495) 
Op 
OM) 24 (D-496) 
oq 
a) =0 (D-497) 
om) _ 1 Z 
OO = 2 (D-498) 
a) =0 (D-499) 
OM) _ ; 
“oa = ° (D-500) 
OM) _ : 
ae oe (D-501) 
O(M) _ : 
= 0 (D-502) 
a(M) _ : 
CD a0 (D-503) 
O(M) x! Vo 0.7p0 oT 7 
dh? Premera. dh Ca 
AM) 4 (D-505) 
Ox 
oO) = (D-506) 
oy 
at) oe 
OG (D-508) 
ag 
a(M) _ 
ag = 8 (D-509) 


D.3.15 Reynolds number derivatives.— 


d(Re) 


3 oe (D-510) 
ae = 0 (D-511) 
oRe) =, (D-512) 
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D.3.16 Reynolds number per 


OV Le 
(Re) _ 
0a et 
O(Re) _ 
oa 
d(Re) 
oe 
O(Re) _ 
a 
O(Re) aii 
Ow 
O(RRe) _ Vat Bp pVol Oy 
Oho” at Oh ow? Oh 
O(Re) _ 
Oz sae 
O(Re) 
ay aca 
O(Re) _ 
oa 
O(Re) 0 
ap 
O(Re) _ 0 
06; 


(D-513) 
(D-514) 
(D-515) 
(D-516) 
(D-517) 
(D-518) 
(D-519) 
(D-520) 
(D-521) 
(D-522) 
(D-523) 


(D-524) 


(D-525) 
(D-526) 
(D-527) 
(D-528) 
(D-529) 
(D-530) 
(D-531) 
(D-532) 


(D-533) 


A(Re’) _ Vo Ip 


Oh ue Oh 
O(Re’) 25h 
Ox 
O(Re’) _ 0 
dy 
8(Re’) 2a 
Oa 
O(Re’) 0 
0p 
O(Re’) = 
06; 


AG) _ 9 
“Op 

AF) _ 4 
“q_ 

A) _ 9 
“Or. 

ag) _ Va 
OV 

AD _ 
Oa 

a) _ 
Op 

AD _ 4 
“Ob 

a@ _ 4 
“00 

AG) _ 4 
He Ve Op 
“Oh 2 
AM) _ 9 
Ox 

Aq) _ y 
“Oy 

AG) _ 4 
0a 

@ _ 9 
OB 

a) _ 4 
05; 


(D-534) 
(D-535) 
(D-536) 
(D-537) 
(D-538) 


(D-539) 


(D-540) 
(D-541) 
(D-542) 
(D-543) 
(D-544) 
(D-545) 
(D-546) 
(D-547) 
(D-548) 
(D-549) 
(D-550) 
(D-551) 
(D-552) 
(D-553) 


(D-554) 
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D.3.18 Impact pressure derivatives.— 


14a M7(1.0 + 0.2M2)25 


Da 5.76M2_\?° 
@ [2-4 (ace va) 


+ 3.0M? ( 


PaVo 
a’ 


5.76M2 ) 9.216M 
5.6M?-08/ (5.6M? — 0.8)? 


5.76M2_\7> 
[pam (528 - 75) 


+ 3.0M? ( 


_5.76M? _ 
5.6M* — 0.8 


. 


(M < 1.0) 


9.216M 
(5.6M2 — 0.8)? 


| 


(M > 1.0) 


Oa 
Oh 


(M > 1.0) 


(D-555) 
(D-556) 


(D-557) 


(D-558) 


(D-559) 
(D-560) 
(D-561) 
(D-562) 


(D-563) 


(D-564) 


(D-565) 
(D-566) 
(D-567) 
(D-568) 


(D-569) 


D.3.19 Mach meter calibration ratio derivatives.— 


Atta oe (D-570) 
acl) =90 _ (D-571) 
Agel) = (D-572) 
14 4(1.0 + 0.2M?)25 (M < 1.0) 
a ae ie (sits) ere 
+300 (SBMS) oaitttae] (2 10 
Atel) ae (D-574) 
Mec aay (D-575) 
ital =0 (D-576) 
O(ae/Ps) =, (D-577) 
Maia) 2g (D-578) 
— 140M (1.0 + 0.2M?)?5 92 (M < 1.0) 
teh Pa) 2 ~ 4 [eam (<a) (D-579) 
+ 3.0M? (ee) Op Ga (M > 1.0) 
Pae/ a) =0 (D-580) 
“acim =, (D-581) 
Aisle) a% (D-582) 
“tsciPe) =H (D-583) 
ae =0 (D-584) 
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ny 


D.3.20 Total temperature derivatives.— 


(Tt) 


= 9 (D-585) 
am) be (D-586) 
ee = 0 (D-587) 
af) : ca (D-588) 
a =% (D-589) 
ve =0 (D-590) 
ae = {10 Sa ae aa oa aa } om eae) 
ah) =0 (D-592) 
a =0 (D-593) 
a) =0 (D-594) 
aa =0 (D-595) 
=0 (D-596) 
D.3.21  Flightpath angle derivatives — 
ne =0 (D-597) 
We =0 (D-598) 
2) =0 (D-599) 
o0) = eae a (D-600) 
ao) _9 (D-601) 
OP 2% (D-602) 
ae aay (D-603) 
OGD) 2g (D-604) 


OG) 2.6 (D-605) 


Ow 
ay) =0 (D-606) 
a) =0 (D-607) 
a =0 , (D-608) 
a) =0 (D-609) 
a =p (D-610) 
a ~ COE (D-611) 
“ = (D-612) 
D.3.22 Flightpath acceleration derivatives — 
ae =0 (D-613) 
ee =0 (D-614) 
ope) =0 (D-615) 
Are =0 | (D-616) 
ee) = 0 (D-617) 
=0 (D-618) 
a =0 (D-619) 
one) =0 (D-620) 
ore) = 0 (D-621) 
as) =H (D-622) 
2s) =p | (D-623) 
ee) =0 - (D-624) 
ae = ; (D-625) 
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othe)" (D-626) 


=0 (D-627) 


=0 (D-628) 


D.3.23 Vertical acceleration derivatives.— 


ach) _ 


Op 


O(h) 
Og 


qb 
2Vogom 


[sin 89(— cos ag Cp, + sin ao CL) 


— sin gg cos Cy, + cos go cos Mp (sin ag Cp, + cos ag CL,)] (D-629) 


Feelsin 8 (— cos a Cp, + sin ao Ci) — sin ¢o cos 8p Cy, 


+ cos go cos M% (sin a9 Cp, + cos ag Cr, )I (D-630) 
qsb ; : 
Veo 9 (— cos ag Crp, + sin ao CL.) — sin do cos 8p Cy, 
+ cos ¢o Cos 9 (sin ao Cp, + cos a C\,)] (D-631) 
1 ' a ‘ _ OX: 
aati sin 00 [-s cos a9 (pVoCp + GCp,,) + Ssin ag (pVOCL + GCL) + | 
— sin do cos Oo (spvecy + qSCy, + i) 
+ cos ¢o cos Ig [s sin ao (pVoCp + GCp,) 
= O27 
+ S cos a (pVoCy + GCL) — ay : (D-632) 
1 F z Lee OX 
peal sin 09 [-25 cos ap (Cp, — CL) + 5 sin ag (CL, + Cp) + al 
: iS OY; 
— sin gp cos I (ascv, + 5c) 
+ cos $g cos Oo [zs sin ao (Cp, — CL) + GS cos ao (CL, + Cp) — I} (D-633) 
1 
ao [sin 90 (-2s cos ao Cp, + GS sin ao Cig + a) — sin dg cos A (ascv, + “2) 
‘ zZ 
+ cos do cos A (as sin a9 Co, + G9 cos ao Cig t+ | (D-634) 
— Gy, COS dp Cos Op + a4 sin do Cos Ay (D-635) 
Azq COS Ay + ayy sin Bg sin do + az, co8 dy sin Oo (D-636) 
0 (D-637) 


: 1 6) a . 1 Op _ 
= oa sin 0 [- COS Qo (Ssvecp P + #SCp,) + sin ag (Gsveeu oh + SC, )| 
: 1 Op. 
—sin do cos Oo (Ssvecy an + a5Cy,) 


+ cos $o cos Oo [sin a9 (5svéep 9p + 15Cp,) 


oh 

+ cos ag (Gsvec. re + 15Ci,)| } (D-638) 
st) = (D-639) 
ae) ah (D-640) 

2) = Speelsin 9 (~~ cos dy Cp, + sin ao CL,) — sin do cos Ip Cy, 
+ cos do cos A (sin ao Cp, + cosa CL; )] (D-641) 

a a elsin A {— cos ao Cp, + sin ao Ci,) — sin do cos Ao Cy, 
+ cos $9 cos Ip (sin ag Cp, + cosa Cr,)I (D-642) 

8) = al sin A [as(- COS Qo Cbs, + sin ao Cr.) + oe — sin dg cos 89 (ascv, + Fe) 
+ cos do cos A [7s¢sia a9 Cros, + cos ao Cry.) = al \ (D-643) 
D.3.24 Specific energy derivatives.— 

o =0 (D-644) 
=0 (D-645) 
oe) =0 (D-646) 
aCe) = 2 (D-647) 
oe =0 (D-648) 
a =0 (D-649) 
ae! = 0 (D-650) 
a =0 | (D-651) 
ree = 0 (D-652) 
ee) = 4 (D-653) 
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D.3.25 Specific power derivatives.— 


(D-654) 
(D-655) 
(D-656) 
(D-657) 


(D-658) 


(D-659) 
(D-660) 
(D-661) 
(D-662) 
(D-663) 
(D-664) 
(D-665) 
(D-666) 
(D-667) 
(D-668) 
(D-669) 
(D-670) 
(D-671) 
(D-672) 
(D-673) 


(D-674) 


O(N) — gSb : 
pach uly See aV, (0° ag CL, + sina Cp,) 


Op 
O(N GSE ; 
am = 27 (°° ao CL, + sin ap Cp,) 
aN) _ gsb 
roe = OV, (cos ao CL, + sin ap Cp,.) 
O(N 
oe = S|cos ao (pVoCL + GCL) + sin ao (pVoCd + GCp,,)] 
ae = g5(cosap Cl, + sin a9 Cp, — sin ag Cl, + cos ag Crp) 
O(N 2 ' 
“ = q5(cosag CL, + sin a9 Cp.) 
aN) _ 

ag 
O(N) _ 

oS 
aN) _ 

Oy 
O(N 1 op | F Op 
oe =$§ [cos a9 (Gvec. Dh + Cin) + sin ao (5 VeCp => Dh + Cp, )| 
O(n) _ 

ox e 
ON) _ 
Oy 
O(N 
aN) = Me GE (ob ag CL, + sinao Cp, ) 
om = T(cos ag Cig + sin a Cb,) 
ao = G5(cos ao CL,, + sin ao Cp,.) 

+ 


D.3.27 Axial force derivatives. — 


O(A GSb 

ae ae (— sin ao CL, + cos a Cp,) 
O(A GSt 

= % (—sin a9 CL, + cos ao Cp,) 
O(A 

4) = Be sin ao CL, + cos ap Cp, ) 


(D-675) 


(D-676) 
(D-677) 
(D-678) 
(D-679) 
(D-680) 
(D-681) 
(D-682) 
(D-683) 
(D-684) 
(D-685) 
(D-686) 
(D-687) 
(D-688) 
(D-689) 


(D-690) 


(D-691) 
(D-692) 


(D-693) 
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= = S[-sin ag (eVoCL + GCL, ) + cos a0 (PVoCD + GCd,)] 


—— = 95(-sinag CL, + cosa Cp, — cos ag CL — sin ag Cp) 


——- = 0 


SEES 


OD 4 


= § [- sin Qo (Svecu ae + Cus) + cos ag (Sve 


= @5(—sin ao CL, + cosao Cog) 


Op 


h Oh 


— = Fe (—sin ag CL, + cos ag Cp,) 
Sb : 
= oe sin ao CL, + cos ao Cp,) 


= @5(—sin ag Crs. + COs ag Cg.) 


D.3.28 x body axis rate derivatives.— 


bo Se ee | 


= COS ag C08 fp 
= —Vosin ao cos Go 


= —Vocosagsin fo 


+ Co, )| 


(D-694) 
(D-695) 
(D-696) 
(D-697) 
(D-698) 
(D-699) 
(D-700) 
(D-701) 
(D-702) 
(D-703) 
(D-704) 


(D-705) 


(D-706) 
(D-707) 
(D-708) 
(D-709) 
(D-710) 
(D-711) 
(D-712) 
(D-713) 


(D-714) 


aE 8 

A) = 0 

xo 0 

ou 

oy 4 

a6 

Ae = 0 

D.3.29  y body axis rate derivatives — 

a 
2) 
ae = Vocos fig 
20) “i 
av) 
00 
21) 4 
an 
20) 
a> =0 
a = 


(D-715) 
(D-716) 
(D-717) 
(D-718) 
(D-719) 


(D-720) 


(D-721) 
(D-722) 
(D-723) 
(D-724) 
(D-725) 
(D-726) 
(D-727) 
(D-728) 
(D-729) 
(D-730) 
(D-731) 
(D-732) 
(D-733) 
(D-734) 


(D-735) 


89 


90 


D.3.30 z body axis rate derivatives.— 


= sin ag cos Bo 
= Vo cos ay cos Bo 


= —Vosin ag sin fo 


D.3.31 x body axis acceleration derivatives.— 


OD) 


g5b 
2Vem 

gGSc 
2Vom 

g5b 


~ 3Vom 


1 


(— cos a9 Cp, + sin ag CL, ) 
(— cos ag Cp, + sin a9 CL,) — Vo sin a cos fo 
(— cos ag Cp, + sin ao CL,) + Vosin Bo 


= {sI- cos ao (pVoC + GCp,) + sin ao (pVoCL + GCLy)] + 


+rosin Bo — go sin ag cos Bo 


(D-736) 
(D-737) 
(D-738) 
(D-739) 
(D-740) 
(D-741) 
(D-742) 
(D-743) 
(D-744) 
(D-745) 
(D-746) 
(D-747) 
(D-748) 
(D-749) 


(D-750) 


(D-751) 
(D-752) 


(D-753) 


(D-754) 


OX: 
= = |as(- cos ag Cp, + sina CL, + sin ao Cp + cos ap Cy) + “| 


— goVo COS Og Cos Bo 


1 


: Ox 
a [75(~ cos a0 Cp, + sin ao Cig) + = 


OB. 


+ ToVo cos Bo + goVo cos ao sin Bo 


—g cos Oo 
0 
= 1v2¢p 2 4 ac ) +sinan (VRC 32 + oCi.)| 
ml ee Be ne OL gol an teh 
0 
0 
a (— cos ap Cp, + sinap CL;) 
_ g5b ‘ 
= 2Vem — 08 a9 Cp, + sin ao Ci) 


os = = [as(- cos a9 Crp,, + sin a CLs.) + ed 


D.3.32 y body axis acceleration derivatives.— 


ae) 


gSb 


= ~——Cy, + Vosin ao cos fo 


2Vom 
qSe 
2Vom CY, 
qSb 


= = Cy, — Vo cos ag cos Bo 


aos 


—_ = = | s(evecv +qCyy,) a8 oF + po Sin a Cos Bo — 79 COS ag COs Bo 


Yr 
= = ~ (a8¢y, + or + poVo Cos ag Cos 39 + ToVo Sin ag Cos Go 


a 


ices ~ (a8Cv, + a) — poVosin ag sin 89 — roVo cos ag sin Bo 


—=— = gcos 6 cos do 


= —gsin 0 sin do 


w+ -=9 


(D-755) 


(D-756) 
(D-757) 


(D-758) 
(D-759) 
(D-760) 
(D-761) 
(D-762) 
(D-763) 
(D-764) 


(D-765) 


(D-766) 
(D-767) 
(D-768) 
(D-769) 
(D-770) 
(D-771) 
(D-772) 
(D-773) 


(D-774) 
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D.3.33 


1 f_ : 
— = > [as(- SIN Qo CDz — COS Qo Che) + 


(6) fly, Op - 
Bho = (5% OY ap + aCy) 
as) _ 

oe 

O(d) _ 

dy ? 

A(b) qe 

aa BVom° %* 

a(b) gb 

dB 2Vom Ya 

OO) 1 Yr 
Bir = aq (95C¥s, + a 


z body axis acceleration derivatives. — 


qSb 


—— — sin dg Cp, — cos ag Cy.) — Vo sin 
3Vem oCp, 0 CL,) — Vosin Bo 


Sz 
—— = ie — sin a Cp, — cos ag CL,) + Va cos ag cos Bo 
0 
GSb 
= sel sin a Cp, — C08 a CL,) 
0 


1 . 7 
aie ae l-s sin a9 (pVoCp + GCpy,) — Scos a (pVoC + GCL) + 


+ Go COS ag cos Bo — po sin Bo 


1 
We Pas [as(- sind Cp, — cos ap CL, — cosag Cp + sin ag CL) + 


— qoVo sin a9 cos Bo 
OZy 


ap 


= —g cos A sin do 


= ~gsin & cos dg 


= Vem sin ag Cp, — cos ao CL, ) 


b : 
Vem sin ag Ch, — cos ag CL,) 


; 1 1 
- sin a9 (Svc e + Cos) — COS Ag (Svc. g 


OLy 


OV 


O2T 


da 


| 


+e) 


| 


— qoVo cos ap sin 39 — po Vo cos Bo 


(D-775) 
(D-776) 
(D-777) 
(D-778) 
(D-779) 


(D-780) 


(D-781) 
(D-782) 


(D-783) 


(D-784) 


(D-785) 
(D-786) 


(D-787) 
(D-788) 
(D-789) 
(D-790) 
(D-791) 
(D-792) 
(D-793) 


(D-794) 


i 1 : OZ 
7 ee [zs sin Qo Cp,, — cos ag CL.) + rie (D-795) 


D.3.34 Angle-of-attack sensor output derivatives — 


ae) ee -# (D-796) 
wee = eo (D-797) 
aed) 2 | (D-798) 
Mes) an a (D-799) 
teu) ag (D-800) 
Magl =0 (D-801) 
a a | (D-802) 
He) =9 (D-803) 
Mou) 4 (D-804) 
“G0 = (D-805) 
Hos) = (D-806) 
ae 24 (D-807) 
Bes) ah (D-808) 
oe =0 (D-809) 
“ce =0 (D-810) 
D.3.35 Angle-of-sideslip sensor output derivatives. — 
ee iene 2 (D-811) 
“ee 2% (D-812) 
“eel Li 2 (D-813) 
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94 


= Tore — PozB 
av V2 


oh = 
0(B i) 

Ox q 
Obi) _ 
ay q 
OB i) _ 
pao 
(Bi) _ 0 
op 
O(Bi) _ 
mae 


= —Yyn cos do cos Og + zp, sin do cos Oo 


(hi) 
06 

Ohya) _ 
oe 


= £7,008 Oo + yn sin do sin Oo + 2), cos do sin Oy 


(D-814) 
(D-815) 
(D-816) 
(D-817) 
(D-818) 
(D-819) 
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